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THEORETICAL PREDICTION OF THICK WING AND 


PYLON- FUSELAGE-FANPOD - NACELLE 
AERODYNAMIC CHARACTERISTICS 
AT SUBCRITICAL SPEEDS 


PART I - THEORY AND RESULTS 


By J.R. Tulinius 
Los Angeles Aircraft Division 
Rockwell International 


SUMMARY 


This report describes the theoretical development and the comparison of 
results with data of a thick wing and pylon-fuselage-fanpod-nacelle analysis. 

The analysis utilizes potential flow theory to compute the surface velocities 
and pressures, section lift and center of pressure, and the total configuration 
lift, moment, and vortex drag. The skin friction drag is also estimated in 
the analysis. 

The perturbation velocities induced by the wing and pylon, fuselage and 
fanpod, and nacelle are represented by source and vortex lattices, quadri- 
lateral vortices, and source frustums, respectively. The strengths of these 
singularities are solved for simultaneously including all interference effects. 

The wing and pylon planforms, twists, cambers, and thickness distributions, 
and the fuselage and fanpod geor 5 tries can be arbitrary in shape, provided 
the surface gradients are smooth. The flow through nacelle is assumed to be 
axi symmetric. An axisymmetric center engine hub can also be included. The 
pylon and nacelle can be attached to the wing, fuselage, or fanpod. 

The wing can have both leading and trailing edge plain flaps of either 
the full or partial span type. Subcritical compressibility is accounted for 
by means of the second order correction developed by Labrujere. 
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INTRODUCTION 


In order to improve the cruise efficiency of v/STOL aircraft with fan- 
pods, pylons, and nacelles it is necessary to locally shape the configuration 
to minimize component induced adverse pressure gradients which cause premature 
shock induced separation. The wing-fanpod juncture must be contoured to 
minimize fanpod on wing interference. Also, the nacelle and pylon must be 
located such as to minimize pylon and nacelle interference on the wing and 
fanpod. If the surface velocities and pressure coefficients can be computed, 
with all interference effects included, then the optimum location and shape 
of the pylon, fanpod, and nacelle can be determined. 

In this analysis the surface velocities and pressure coefficients are 
computed for a thick wing and pylon-fuselage -fanpod -nacelle combination. The 
wing and pylon induced flow fields are represented by source and vortex 
lattices. The source lattice induces flow due to airfoil thickness and the 
vortex lattice induces flow due to lift. The flow induced by the wing and 
pylon thickness can be computed directly since the source lattices are placed 
on the chordal planes, which allows the source strengths to be defined by the 
local slopes of the distribution of wing and pylon thickness. The wing and 
pylon vortex strengths, however, must be determined simultaneously along 
with the quadrilateral vortex strengths for the fanpod and fuselage and the 
source frustum strengths for the nacelle. 

In order to determine the wing, pylon, fuselage, and fanpod vortex 
strengths and the nacelle source frustum strengths a system of aerodynamic 
influence equations is developed and solved. The equations are developed by 
computing the amount of velocity induced by each of the vortices and source 
frustums, of unit strength, normal to the surface of each component at a 
series of control points. 'Tiis set of influence equations is reduced in size 
by constraining the strengths of the vortices in the chordwise direction by a 
vorticity series. The equations are then solved by partitioning the matrix 
into two sets of equations: l) that due to the vortex strengths, which is 

overdete mined, and therefore solved by a least squares technique and 2) that 
due to the source frustums. Both inversions in the solution of these two 
sets of simultaneous equations are done by Householder's procedure. The 
program has options which allow the smaller interference influence matrices 
to be neglected in the solution, in order to save computer time. If these 
influence matrices are to be neglected they are neither computed nor manipu- 
lated in the solution of the influence equations. 

In this analysis advantage is taken of the similarities between the 
source and vortex induced flow fields. The influence equations for the source 
and vortex lattices have been formulated in terms of the same expressions. 

This permits the velocity induced by the source lattice and that induced by 
the vortex lattice to be computed and summed simultaneously, saving approxi- 
mately 50 percent of the computing time that would otherwise be required to 
calculate the wing and pylon influence matrices. Also, since quadrilateral 
vortices are used on the fanpod and fuselage, it was possible to neglect the 
contribution of those vortices which are more than a given distance from a 
control point. The induced flow from a quadrilateral vortex decreases rapidly 
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as the control point is moved away from the vortex. This is due to the 
cancelling effect of adjacent parallel sides. The induced flow decreases 
proportional to (~)^ where A is the perpendicular distance between two 
adjacent parallel sides and T is the distance between the centroid of the 
vortex and the control point. A similar effect is obtained in the spanwise 
direction on the wing and pylon due to the cancelling effect of adjacent 
trailing ortex legs. Because of these effects } limits were put into the 
program which define the significant area of influence of a vortex element and 
contributions are computed only for those control points which fall within 
this area of influence. The limits are given as input data. This concept has 
produced a considerable savings of computer time with no significant change 
in accuracy. 



LIST OF SYMBOLS 


A consistent set of units is assumed throughout this report excepu for 
the skin friction drag equations derived in Appendix K, which assumes the 
units specified in that appendix. 


A 

A 

b 

f 

C 

Cd 

CD 

cl 

Cn 

Cm 

cm 

Cp 

a 

AS 
V S 

E 

F 

h 


Aerodynamic influence matrix 
In remental area 
Angle of attack 
Span 

Pylon angle of cant or compressibility factor 
Chord 

Section drag coefficient 
Total drag coefficient 
Section lift coefficient 
Total lift coefficient 
Section normal force coefficient 
Section pitching moment coefficient 
Total pitching moment coefficient 
Pressure coefficient 

Vorticity, ratio of specific heats, or pylon dihedral angle 

Fanpod, fuselage, or root section vortex grid length or width 

Freestream plus source lattice flow normal to component 
surface 

Influence function 

Wing or pylon chordwise vorticity series coefficient 
Pylon height 
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Fanpod or fuselage polar angle measured from the X-Z plane 
about the X axis. The positive direction is clockwise when 
looking from the fanpod tail end to the nose. 

Discrete vortex strength 

Free stream Mach number 

Unit vector normal to surface 

Number of vortices and sources in chordvise direction on the 
wing or pylon 

Number of terms in the chordvise vorticity series on the 
wing or pylon 

2^w/bw or Y P/h 

Special spanvise load shape or J/ / 

Chordvise polar coordinate equal to COS "^1-Cx/c^) or sweep 
angle 

Influence matrix due tc source lattice 

Surface source strength 

Unit vector tangent to surface 

Perturbation velocity in X direction 

Total velocity 

Free stream velocity 

Perturbation velocity in Y direction 

Perturbation velocity in Z dl. 'ction or width 

Spanwise source lattice spacing 

Spanvise vortex lattice spacing 

Camber line ordinate 


Wing or pylon thickness 


AR 

S' 
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A 


u/ 


i vvj 


wur 


^Ttvir 


•tv. 




U 




u. 


ur^ 


•*y 


Aspect ratio b 2 /s w 
Flap deflection 
Wing twist 

Local sweepback angle of constant $ chordlines 

Velocity components in equations (129) through (138) 

Total surface perturbation velocity in chorduise direction 
including lift and thickness interaction 

Total surface perturbation velocity in spanwise direction 
including lift and thickness interactions 

Chordvise component of perturbation velocity due to lift 

Spanwise component of perturbation Velocity due to lift 

Chordwise component of perturbation velocity due to thickness 

Additional thickness perturbation velocity term which does 
influence lift 

Spanwise component of perturbation velocity due to thickness 

Additional thickness perturbation velocity term which does 
influence lift 
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Subscript 

B 

BB 

BF 

BIT 

BP 

BW 

C 

F 

f 

FB 

FF 

FN 

FP 

FIT 

i 

tso 

J 

NF 

MN 

HP 

nw 

PF 

PW 

PF 


Definition 

Fuse luge 

Fuselage on fuselage 
Fanpod on fuselage 
Nacelle on fuselage 
Pylon on fuselage 
Wing on fuselage 
Chord or camber 
Fanpod 
Flap 

Fuselage on fanpod 
Fanpod on fanpod 
Nacelle on fanpod 
Pylon on fanpod 
Wing on fanpod 

Source or vortex location in chord-wise direction, inboard 
lattice panel, or incompressible coordinate. 

Freestream direction 

Control point location 

Fanpod on nacelle 

Nacelle on nacelle 

Pylon on nacelle 

Wing on nacelle 

Fanpod on pylon 

Nacelle on pylon 

Pylon on pylon 
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J 


PW 

WF 

WN 

WP 

ww 

K 

L 


L.E. 


M 

N 

0 

P 

T 

t 

u 

V 

V 

(XN,Yn) 

(Xp ? Yp,Zp) 

(Xw,Yy,Z w ) 

(Xb,Y b ,Z b ) 


Wing on pylon 
Fanpod on wing 
Nacelle on wing 
Pylon on wing 
Wing on wing 

Leading edge flap or span location index 
Spanwise location of source or vortex 
Leading edge 
Longitudinal direction 

Nacelle, normal direction, chordwise load shape 

Outboard lattice panel 

Pylon 

Lateral direction 
Thickness 

X direction perturbation velocity 
Y direction perturbation velocity 
Z direction perturbation velocity 
Fanpod coordinate system 
Nacelle coordinate system 
Pylon coordinate system 
Wing coordinate system 
Fuselage coordinate system 
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THEORETICAL DEVELOPMENT 


CONFIGURATION REPRESENTATION 

The theory discussed in this report is capable of predicting surface 
velocities and static pressure coefficients, section loads, and total loads 
and moments for a configuration cf the type shown in figure 1. 




(0,0,0) 


Figure 1. Thick Wing and Pyl on- Fusel age - 
Fanpod- Nacelle Configuration 


In general, the fanpod and fuselage can he arbitrary in shape. The 
fanpod can be placed anywhere on the wing except at the wing tip. The fanpod 
and fuselage do not have to have straight juncture lines where the wing and 
pylon are attached. This allows the fanpod and fuselage to be shaped in any 
desired manner, in the juncture region, to produce a favorable induced 
effect on the wing and pylon. In order to do this the wing and pylon have 
been divided into two sections; the basic wing or pylon panel and the root 
section. The basic wing or pylon panel is assumed to be near planar. 
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The root section is a transition region from the basic panel to the juncture 
line where the wing or pylon is attached. 

The wing and pylon can be arbitrary in planform shape and can have any 
desired twist, camber, and thickness distribution. The wing can have either 
full or partial span leading or trailing edge flaps. 

The nacelle is assumed to be axisymmetric and is attached to the pylon. 

A center engine hub can also be included with the nacelle. Also, different 
mass flow ratios can be simulated by blocking different amounts of flow 
going through the nacelle. 

The perturbation velocity induced by the fanpod and fuselage is produced 
by a grid of quadrilateral vortices on the external surface of the fanpod 
and fuselage. A quadrilateral vortex is a line vortex of constant strength 
bert into the shape of a quadrilateral. As shown in figure 2 vortex lines 
are placed both collinear with and perpendicular to meridian lines. The 
strengths of these lines are assumed to be constant around each of the quadri- 
laterals formed by the intersection of longitudinal and lateral grid lines. 

The strength of any one vortex segment is equal to the sum of the strengths 
of the two adjacent quadrilateral vortices. The direction of the circulation 
due to a vortex segment is defined by the right hand rule when applied about 
the segment. The positive direction is assumed to be clockwise around the 
perimeter of the quadrilateral when looking at its external surface. 



Figure 2. Fanpod Quadrilateral Vortex Grid 


lU 



Figure 3. Fuselage Quadrilateral Vortex Grid 


The flow induced by the basic panels of the wing and pylon is produced 
by source and vortex lattices placed on the wing and pylon chordal planes. 

As shorn in figures 3 and ^ the bound vortex and source lines are skewed 
and collinear with each other. The trailing vortex lines are in the free- 
stream direction. The bound vortex and source lines are spaced at equal 

increments of or <sv , where J and 

C£>s‘l> -z<*r- , This type of spacing places more singularities 

near the leading edge where the source and vortex strengths vary most rapidly. 
The trailing vortex lines are evenly spaced. 



Figure it. Wing Source-Vortex Lattice 
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The flow induced by the wing and pylon root sections is produced by 
quadrilateral vortices and skewed source lines. As with the basic panel, 
the vortex and source lines produce flows which are antisymmetric and sym- 
metric about the mean camber surface, respectively. The antisymmetric flow 
represents the solution due to lift and the symmetric flow represents that 
due to thickness. The quadrilateral vortices in the root section are formed 
by the intersection of lines equally spaced between the juncture line and the 
basic panel, and lines spaced at equal increments of u or in the chord- 
wise direction. The skewed source lines are collinear with the bound vortex 
lines as in the basic panel. 



Figure 5. Pylon Source -Vortex Lattice 


The pylon vortex lattice extends up into the nacelle in order to satisfy 
the boundary condition of no cross flow through the nacelle. This is nec- 
essary since the nacelle is assumed to be axisymmetrie, which makes it 
possible to satisfy the nacelle boundary conditions at only one lateral 
station. These are satisfied at the lateral station closest to the surface 
where the pylon is attached, to ensure the best possible representation of 
the nacelle interference within the axisymmetrie nacelle assumption. 
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Figure 6. Flaw Through ife.ce.lle Source Frurtum Representation 

'The nacelle induced flow is produced by a distribution of source frustums 
as shown in figure 6. The frustums are distributed around the inner and 
outer surfaces of a flow through nacelle or over the external surface of any- 
desired axisymmetric shape. A center engine hub or other body used to block 
mass flow can also be represented by source frustums along with the nacelle. 


INFLUENCE EQUATIONS 

Each of the skewed sources, source frustums, quadrilateral vortices, and 
skewed vortices induce flow at control points distributed over the surface of 
the thick wing and pylon-fuselage-fanpod-nacelle configuration. Their 
strengths are solved for simultaneously such that the flow induced by them 
plus the free stream, is a minimum normal to the configuration surface at the 
control points. In order to obtain the singularity strengths a set of aero- 
dynamic influence equations is developed by computing the flow induced normal 
to the surface at each control point due to a unit strength of each singularity. 
These are then multiplied by their respective unknown coefficients, summed 
along with the component of frees tream velocity normal to the surface, and 
set equal to zero at each control point. This results in a system of equations 
where the number of unknowns are equal to the number of singularities and 
the number of equations are equal to the number of control points. 

The number of unknowns can be reduced by constraining the singularity 
.'brengths by a set of finite series. When this is done, the number of unknowns 
cc comes equal to the number of terms in the series times the number of series 
used. In this analysis the skewed vortices on the wing and pylon and the 
quadrilateral vortices on the fuselage and the fanpod are constrained in the 
chordwise direction. Also, there is the option of constraining the skewed 
vortices in the spanwise direction on the wing. When these singularity 
strengths are constrained the aerodynamic influence equations described above 
are linearly transformed such that the unknowns are then the coefficients 
of the terms in the finite series. This results in a system of equations 
which is overdetermined. The equations are solved using a least squares 
technique. 
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WING ON WING 


The influence of the wing on itself is divided into two parts; the 
influence of the vortices and the influence of the sources. The sources are 
assumed to lie on a plane and can therefore only produce velocity normal to 
that plane at their respective locations. Because of this fact, the "wing 
sources cannot affect the strength of each other nor can they affect the 
strength of the wing vortices. Their influence can therefore "be omitted 
from the wing on wing influence equations. 

The velocity induced hy the wing vortices normal to the wing mean camber 
surface is given by; 

{^7} = [**. ^ } &) 

Where and /\£, 3 are defined in the following sketch. 



The local angle is assumed small and the wing dihedral zeroj so 
that j m* =« soi ^ = o ana aV. / 

Therefore ; 

=■ {[>%][^] + ^ 

Since; both and are small the product P' 

will he omitted. 

Therefore ; 

C3) 

The elements of /&'*''**, axe computed with equations for the velocity 
induced hy skewed and quadrilateral vortices } defined in Appendices C and D, 
respectively. 
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FANPOD ON WING 


The fanpod is represented by a grid of quadrilateral vortices on the 
fanpod surface. The velocity induced by the fanpod quadrilateral vortices 
normal to the Tiring mean camber surface is given by 


{t'}“ [ 


fe} 

(4) 

The influence matrices /&*, 

, /?««■ , and /?«,- are defined as: 

V rw 



i * 


J 

(5) 

it 

1 J 

f j 



(6) 

l A ~"*~ ] =-| 



(7) 

Where is defined in the 

following sketch. 





The angle and the matrix are assumed small, so that 

cosg-x-r^ ■sz. / j and the product ] 

is omitted. 

Therefore : 

fe} ‘ {['**" J 

+ [ Jjfej (8) 
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/ 



In addition, since , and are small, the products 

WH1 and j are omitted. 


Therefore ; 


{%} = ["Kit} 


( 9 ) 


The elements of are computed with equations defined in Appendix D. 

«J= 


FUSELAGE ON TONG 

The influence equation for the fuselage on the wing is analogous to 
that for the fanpod on the wing. Therefore; 

Pt*} = [ Arn> ** Ml 

The elements of are computed with equations defined in 

Appendix D. 


PYLON ON TONG 

The influence of the pylon on the wing is due to the vortex and source 
lattices on the pylons. The velocity induced by these lattices normal to 
the wing mean camber surface is given byj 

* 'K'K, * K, 4%, W 









( 11 ) 


Where the matrices , 

influences and the matrices 
source influences. 


, and are due to the vortex 






and Sivp- are due to the 
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The pylon can be rotated by an angle Y and positioned at any point 
as shown in the following sketch. 



and 

The influence mattrices 
are given by: 
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Since , and are small, the products 

and f' J s^j 1 are omitted. Also, since the wing has no 
dihedral, ^ =c? " p 

‘f3 

Therefore ; 



The elements of , ^u7j, p , and are computed 

with equations defined in Appendices B, C, and D. The port pylon influence 
is computed "by moving the control point to its image location on the port 
side, computing the influence of the starboard pylon at that point, and then 
changing the sign of the "Y" component. 


NACELLE ON WING 

The flow through nacelle is represented by a system of source frustums 
as shown in figure 5 - The velocity induced by the nacelle frustums normal 


to the wing mean camber surface is given by; 

$$ = + '’f* 'Mw * 7£§sj (19) 

The influence matrices , and are defined as; 

]/W J'fa'tJ (20) 

= " ( 21 ) 

P**J "" ‘ (22) 

©w/j ' d°s . <k J ( 23 ) 
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Again since £ ana j/kwvj are small, their product is omitted 


Therefore ; 


( 25 ) 


The elements of are commuted with equations derived in Appendix E. 

The influence of the port nacelle is computed by moving the control point 
to its image location on the port side, computing the influence of the 
starboard nacelle on that point, and then changing the sign of the "Y" 
component. 


FREESTKEA.M ON WING 


I - f A-L 4- /^by -/■ 1 

luTi L \& *■ or 

where v ^ 

_ o . ✓sWz? — o 

J ~vT J ^ 

Therefore; 

{%} = f ^ J ■ } ■ H 

WING ON PYLON 


(26) 


( 27 ) 


The influence of the wing on the pylon is due to the source and vortex 
lattices on the wing. The velocity induced by these lattices normal to the 
pylon mean camber surface is given by the following expression. 


{%} = ['%>'*** 


+ 


** /S Ns- ]{% } 


( 28 ) 
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Where and s4> are defined in the following sketch. 



If the local angle PPr is assumed small and the pylon cant zero, 
then Ai> gpa j =o , and /£> = to . 

fjc? cjd 

= { ]E ^ ^p] + [ {xit} 


+ 


{wiw -[WM 


(29) 


Since f FEtj , , and jW»*J are small, the products 

are omitted. 

Then considering the difference in orientation between the wing and 
pylon coordinate frames : 


[*v] 


( 30 ) 

[W 


( 31 ) 


Therefore ; 




+ wjwiife} 


(32) 


The elements of the matrices , and 

are computed with equations defined in Appendices B, C, and D. 
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FANPOD ON PYLON 


The velocity induced "by the fanpod quadrilateral vortices normal to 
the pylon mean camber surface is given by: 

{%} * \*K 'H . " ^ 'H. ^ ^ It } 

After taking into account the difference in orientation of the pylon 
and fanpod 


(33) 


['H] 

+ [ Vwjr j[ 


(34) 


(35) 


(36) 


Therefore, since 


4 - J] [t] (37) 

Since both and are small, their product is omitted. 

x f 

(%}"{■ i + r os) 

The elements of and are computed with equations defined 


in Appendix D. 


FUSELAGE ON PYLON 

The influence equation for the fuselage on the pylon is analogous to 
that for the fanpod on the pylon. Therefore; 

rs " [-[ W J[^1 + MMM3 ( 39 ) 
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The elements of ^pb v and. / %s^ are computed with equations defined 
m Appendix D. 


PYLON ON PYLON 

The Influence of the pylon on Itself is computed in the same manner 
as the influence of the wing on itself except for the influence of the 
port pylon on the starboard pylon. 

]{■ ■f'j 

- L'*- **s- - '•K If! (to) 

Both the source and vortex lattices of the port pylon contribute to 
the starboard pylon. Only the vortex lattice of the starboard pylon 
contributes to the starboard pylon influence equations. 

oince; /,1 ^r =c> f ^^^/.(sand , J ~ j and F»jf 

are small 

[%} ' [~ r—JIn 1 + 

+ ]}[^} * 

+ ** I) fe] 

(*a) 

where » and are influences of the port 

pylon onto^the starboard pylon. 

The five level subscript notation in equations (^2) through (^5) 
indicates: Levels 1 and 2) FPp influence of port pylon on starboard pylon; 

Levels 3, k f and 5 ) Xpp coordinate of vortex segment located on port pylon. 
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L’Sl 

= I^j 

1 + 

(42) 

i*%] ■= 

i - 

(43) 

C^j 

I = - 


(44) 

and 


(45) 


The elements of } > and are 

computed wiuh equations defined inftppendices B, C, and D. The influence of 
the port pylon is obtained by moving the control point to its image location 
on the port side and then computing the influence of the starboard pylon 
at that point. When this procedure is used, the sign on the "Y" component 
of velocity must be changed. This sign change is applied in equations (42) 
and (44). 


NACELLE ON PYLON 

The expression for the velocity induced by the nacelle normal to the 
pylon mean camber surface is given by the following expression 

[£&} = [ ^ '^■v ” ^Hs. ][-f } ( 


Since the nacelle produces an axisymmetric flow field, no contribution 
to is obtained f ; om the starboard nacelle on the starboard pylon. 

However, both nacelles contribute to and , Since and 

are small and since M> - o ; 

M = [ 'H. ]£§} cw 
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Where 





and 


L^vl '- 





m 


(<©) 

( 50 ) 


The elements of •^^^are computed from equations derived in Appendix E. 
The influence of the nacelle on the port side is determined by moving the 
control point to the port side and computing the influence of the starboard 
nacelle at that point and then changing the sign on the "Y" component. 


FREESTREAM ON PYLON 
Where 

* 

Ak ,*%=><=> . j ISSP ^ '• ° 

and -=c 

Therefore ; 

WING ON FANPOD 


( 51 ) 


( 52 ) 


The wing influences the fanpod quadrilateral vortex strengths with the 
flow induced by the wing source and vortex lattices. The component of this 
flow normal to the fanpod surface is given by the following expression. 



y- 



■f 


L 


/lA 


■vw. 





r r i^’2. 



( 53 ) 
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where 


I>v] 


JVwpa^j^J 

(54) 

L /?w * ] 

’ Lf^*] 


(55) 

3 


j^z ^ J 

(56) 

[>"», J 



(57) 

i — i 

V 

i 

if 

£ 1 

= [ ^ ] 


(58) 

[ s «v3 



(59) 


Since > and are small; 

{%} * (IHKI - Wsl * KW + 

+ (60) 

1 

Also, since the product i s sroall; 

[%} = [C<M + KM J 

+ (lH 'JL S ^ J * F*'* t Sl) 


The elements of /7,=w XfJ , , and and of , 

and are computed with equations derived in Appendices C, B, and D. 

€ij 

The components of the unit vector normal to the fanpod surface are 
computed from the coordinates of the vortex grid. 
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(62) 


~~ 'Ce_ __ • 7 ' -7— 


/ ^V w ^r - — / -_ 7 — 

^ r Xfl 7^ 


A" "■%=- • 


and. 


where 

^ ~ J h * *** - **,)/as r „. 


^ = ?(% - % +Yr t - Y*,)/a , 




C r* 


= -• " ^ 


">7 • 


%" C ** ~ x "‘ -^J/as^.. 

S.- =? fe - », ' X, - * s j/^ 

/ 'i 


and 


% = ?<r ^> - ^ - 1 ^ - *s y.^. 

Also where; 

^ '**• * ***-**,? *■ v\ \ - y f , + y Pi - v^y- 

■^H 7 C^-? Fi *- - ?PB f j *i 


(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


(71) 
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anu. 



+ 




+ ® v -s* 1 )] f '' t ( 7S ) 


Where the indices refer to the corners of the quadrilateral vortex 
as shown in the following sketch. 



The i subscript indicates that the unit vectors are for the equivalent 
incompressible fanpod obtained from, the Gothert compressibility transformation. 


FANPOD OK FANPOD 

The influence of the fanpod on itself is due to the velocity induced 
by the fanpod quadrilateral vortices normal to the fanpod surface. This 
component of velocity is given by; 

{%} “ fa** J^} 

The elements of , and are computed from 

equations defined in Appendix D. 
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FUSELAGE ON FANPOD 


The influence of the fuselage on the fanpod is analogous to the 
influence of the fanpod on itself* 

Therefore ; 

{^r} ~ ^ ^ ^ ][~%] (7*0 

Where the elements of , and ''Sks^ are computed with 

equations defined in Appendix D. 


PYLON ON FANPOD 


The influence of the pylon on the fanpod is due to the pylon source 
and vortex lattices. These lattices induce a velocity normal to the fanpod 
surface which is given "by; 



When the difference in pylon and fanpod coordinate frame orientation 
is taken into account, 

+ fc*r , J [ 

■* ]j)jv] (76) 
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The elements of j <s: ^% r , and 

are computed with equations defined in Appendices B, C, and D, The influence 
of the port pylon is computed by moving the control point to its image 
location on the port side, computing the influence of the starboard pylon 
at that point, and then changing the sign on the "Y" component. 


NACELLE ON FANFOD 


The velocity induced 
frustums is given by: 


normal to the fanpod surface by the nacelle source 


The influence matrices , and are defined as 

l>*„ J ~ ] 

]_?***& ] ” [ X ^ 


( 77 ) 


( 78 ) 

( 79 ) 
( 60 ) 


where 
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The influence matrices and are computed with equations 

defined in Appendix E. The influence of the port nacelle is computed by 
moving the control point to its image location on the port side, computing 
the influence of the starboard nacelle at that point, and then changing the 
sign of the "Y" component. 


FREESTREAM OI'J FANPOD 

[^t] = [ ^ 

Therefore; 

s “* w, vi{ 1 '} Jf-) m 


WING ON FUSELAGE 

The influences of the different components on the fuselage are computed 
in a manner analogous to that for the fanpod. Therefore; 



The elements of , and 

are computed with equations derived in Appendices C, B, and D. 

The surface unit vectors O’fcg » ^ys ^ ,) 

and (”2s 7 ^ , /sry e are computed the same as those for^the fanpod 

and defined by equations (47) through (55). 


FANPOD ON FUSELAGE 

[ %-] = [ ^ ^ ^ ^ ^ ^ ^ J-f j (86) 

The elements of , and are computed from equations 

defined in Appendix D, 
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FUSELAGE ON FUSELAGE 


]Rfj ( S7 ) 

The elements of , ^3a, , end ''' 3 ^^ are computed from equations 

defined in Appendix D. 

PYLON OM FUSELAGE 

[%] = JKJ + r^fML'Vl 

* [r^i^ ] (e8) 

The elements of > and ^%=> 

are computed with equations defined ir. Appendices B, C, and T3. The influence 
of the port pylon is computed by moving the control point to its image 
location on the port side, computing the influence of the starboard pylon 
at that point, and then changing the sign or. the "Y" component. 


3? 



NACELLE or; FUSELAGE 


{%] 

^KJMK]}[f] (89) 

where 

... s g<>s * - 7 

e s»’ bOt}-y^)-s-.~jf -fy? 4 )s»s*' J (90) 

The influence matrices /%-yU^ and are computed with equations 

defined in Apnendix E. The influence of the port nacelle is computed as in 
the case of the nacelle on fanpod. 

FREESTFEM4 Oil FUSELAGE 

f%} = f [ 1 ] + [ s ^ 1 ] (91) 


WING ON NACELLE 

The flow induced normal tn the nacelle surface by the wing source and 
vortex lattices is given by the following equation: 

+ - 5 ,^ + sr^ ][& } 


After transforming the influence matrices to the wing coordinate system: 
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(93) 



The influence matrices / ‘^Wx uy > ? and are computed with 

equations defined in Appendices C and D, and the matrices , 

and are cormuted with equations defined in Appendix B. Also, 

y.j _ dyj. and sO -= c.os e*t.. where at is defined in Appendix E . 

' <rT V L) ' t 'tjj 


PYLON ON NACELLE 

f Phe flow induced normal to the nacelle surface by the pylon source and 
vortex lattices is given by the following equations. 



After transforming the influence matrices to the pylon coordinate system: 

P8- - r<i[~v ]] [p 


where 


L 


'H, 


1 


+ [ N e oS ^ [ + 

4 ] 


(95) 


(96) 
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ch.j ~ r w v] {■ + ml - a 



+r 

L. ■£=> J 


The influence matrices 3 y , and 5^g.^are due to 

the port nylon and and <r^ are due to the tarboard pylon. The 

influence of the port pylon is computed by moving the control point to its 
image location on the port side, computing the influence at that point, and 
then changing the sign on the "Y" component. 

The influence matrices and are computed from equations 

defined in Appendices C and D, and the matrices and are 

computed vith equations defined in Appendix B. 


FANPOD ON NACELLE 

The flow induced by the fanpod quadrilateral vortices normal to the 
nacelle surface is given by: 

{ ^ ( 98 ) 

After transforming the influence matrices . and / ^%j to the fanpod 
coordinate frame; 



The influence matrices ^4=^ , j and are computed with 

equations defined in Appendix D. 
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FUSELAGE ON NACELLE 






( 100 ) 


The influence matrices 
equations defined in Appendix D. ' 3 


and are computed with 


NACELLE ON NACELLE 

The flow induced by the nacelle frustums normal to the nacelle surface 
is given by: 




The second subs crip 4 : a nd Aj$ refer to the port and starboard 

nacelles , respectively. The port nacelle influence is computed by movinr 
thp control point to its image location on the port side, computing the 
influence at that point, and then changing the sign of the ”Y" component. 

The matrices and t£ksJ Y , are computed with equations defined 

xn Appendix E. 


FREESTREAM ON NACET.LE 



Combined influence equations 


V£© 
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o 

VjpB 

Y V^pp Y 
Ve 
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o 

VVrrjS 
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j. 

VU 
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Vo 
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Vo 

VpJyJ/T) 

Vo 

— o 

V, 

-f- ^PF 
Vo 

-(- Vjpyn) 

V* 

+ yjpp 

Vo 

+ y^'O 

Vo 

•= C 

V^S 

lio ' 
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Vo 

y ^ 

Y 

Vo 

Y~ 

Vo 

■= O 


Therefore, the discrete influence equation is; 



A-er 

/Vv-0 

Asp 

A-p/o 

% 


^3 

Aps 

ApF 

Apyi 

Apr 

App) 

% 


'7‘Sp 

AvS 

Ayjp 

AyJrJ 

ArJT= 

ArirO 

% 

— ■ 

FS W 

/Vs 

A-^p 

AprJ 

/Vf* 

A^pp) 

% 


75p 

A/o& 

A^p 


Amp> 

1 

1 ^ 

% 




(ioh) 

( 105 ) 


(106) 
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CONSTRAINT EQUATIONS 


The discrete matrix influence equation can he transformed to reduce the 
number of unknowns by letting the discrete vortex strengths on the wing and 
pylon be defined by: 


\ [( /<oL /v ao )j' + (*&/&)£ 7 - 

/• J (107) 

A^/ 

Where % , /> > F < , and and l4c , i<$ t - , KU , and /£' are defined in 

Appendices p and G, respectively. Equation 107 can be written in matrix 
notation as: 


fe] a [ ^ p} <“ 

Also, due to the definition of in equation 51 of Appendix F. 

fe} = f T ^) ][ T ^> ]{-} = [ s ]{*} (1C 

The pylon vortex strengths are constrained in the chordwise direction 
only, therefore, they are defined by an expression of the form; 



(no) 


The form of the expression used to define the K's on the wing can be 
either that of equation 108 or that of equation 109, depending on 
whether just the chordwise variation of vortex strengths is constrained or 
both the chordwise and spanvise variations of vortex strengths are constrained. 
If just the chordwise direction is constrained, the wing vortex strengths are 
defined by; 


feM T ~ JH 


(111) 


If both the chordwise and spanvise directions of the wing are constrained; 



( 112 ) 



A matrix constraint equation for the fanpod and fuselage vortex strengths 
Is derived in Appendix H. This equation has the form; 

(U3) 

After combining equations I 10 , 112 , and 113 , the transformation 

matrix by which the discrete influence matrix is post multiplied, to obtain 
the constrained influence matrix, is given by; 

[ 0 ]I o l[°l[ = ] 

!H>] L°H o 3 L°1 (u 'o 

[°n°] i<i l ° ] l ® i 


Therefore, the constrained influence matrix is; 


Asb A^p Aw Asp 


A&s A&r Aw A-gp Asp 

'&L&S & e> CD o 

A\F8 A?? Apr-i App ApsO 


Apg App /}pp -ArpAprJ 

o o o ° 

[ Apa Avjp Aprj AJyjp Aprj 

— 

Apb App AppJ AtppApj/O 

o o o a 

\Afb App ApP> Avp/Uj 


Aps AAp App ApioAfap 

o O o 7p O 

\ App A)pp App App 


Ap& App App App App 

c? <d a & t 
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SOLUTION OF INFLUENCE EQUATIONS 


To solve the constrained influence equation the unknowns are divided 
into two setsj (l) those which define the strengths of the singularities on the 
wing, pylon, fuselage and fanpod, and (2) those which define the source strength 
of the nacelle. When this Is done, the influence matrix equation is parti- 
tioned as follows; 



AbP AK-bvj A&P 

1 


'§■ 


“ — 
^8 

Af& 

A PF Afvi) Afr> jAp/J 


A 


W£T f 

A MB 

ArJF Awia/ A'rJ-p 1 //wrtl 





Ape 

J 

App Afhi ft-p-p | Apri 




VSp_ 


A/op A/o/o Ap/O 


% 




(116) 


This matrix can then he represented hy 


A/ ^/-z 







<ZZ-£ 



— ^ 


_ 


L. — 


where 




'V 


nJ 



/9prJ 







'■'tv'W 


( 117 ) 


(118) 


( 119 ) 





U4 


( 120 ) 



and 


( 121 ) 

The solution to equation 117 is then found as follows; 


/?„ <2, r* ^ (122) 

S9 Z , < 2 / ' ( 123 ) 

^ZZ &-ZJ#; (7:5, (i 2i| ) 

^2^23 * S?sz ^z. = ( 125 ) 


Then after subtracting equation 125 from equation 122; 

(126) 

Since p?„ ~~ J is not square, equation 126 is solved using 
Householder’s triangulation procedure discussed in Appendix A. This solution 
will he denoted by; 

^ R z C , )] [ vs, -s?,* C=s; ) ] (12T) 

Equation 127 is then substituted back into equation 123 to obtain ct e . 

ez* «/■?«' e>^ 


J+5 


(128) 



SURFACE VELOCITIES AMD PRESSURES 


In order to compute the surface velocities, the component influence 
matrices are multiplied “by the components of the surface tangent unit vector, 
in a scalar product, and then in turn post multiplied by the unknown coeffi- 
cients obtained In the simultaneous influence equation solution. This is done 
for the influence of the wing, pylon, fuselage, fanpod, and nacelle on the 
wing, pylon, fuselage, and fanpod. The components of the surface resultant 
velocity due to each of the configuration components are then summed to obtain 
the total surface velocity at each of the control points. This velocity is 
then substituted into the Bernoulli equation to obtain the surface static 
pressures. 

WING OM WING 

The velocity ratio components on the wing due to the wing vorticity and 
source density are given by the following expressions. 






(129) 


0 + o + ± 


'/« 




t ( { J + ^ ]} i fjk 


(130) 


Where the upper sign is applied to the upper surface calculation and 
the lower sign to the lower surface, also; 


(itl-e = <■ a “'*t 


(131) 


(wo ,^‘3" j fJe y ** C&^J i 


+ /*"**<<*> H 

<$,«) lJ 


(132) 


^ V^, locc" z [? ^ + ©U.) - a &7-&/-Z [rs 




(133) 
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and 


. x v 1 r s‘"*o*'-oau 

“ ?C4^,^ ~ ^ f Z_ XT' L ~Z&-,J 

2. ^ ^ 




_ SssiSC*f+Ot s>„ 


7 +. ^_§. f C <s> v 

J *>te„ L 


- ao3<9, \ ^ag /f W ^ g ^gO / (^“S^cS/ 

^ [cacs^ - j + 

w r^'s^^C&n +<&&.) t Jj 


M^r3 

G F (_ / y- ^ £0^ 

" dVS./ 


V °<f’ 


’/ svw^ £< 5 ^ - J / 


*■£ zg><s nf tM g 

fe"*i c^-^/jl 


i(t*^ ?* re - ts?so0£^ X / ~ *«**,)} 


(13 1 *) 


where 

,-f 

^ '*t 

> 15. , and 

■=*'>? 

L u 

Also; 



/«C 

__/ r r 

t4 

***** t ”7 


and z>w *z£ 


and the derivatives 


J 




are defined in Appendix J, 


-<&*{££«. [Psps -ps?-; mi i 


•’it. <• l- "“it. 


~r 

* K> 


*CU‘ . 


(135) 
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- E5&7 { 1 


A."/ £ -f 


J»*u *• G* 

^7. 

Clu 


* . . luJ • 

^HL *- 1 U 


7** 


(136) 




l£o 


/ 


„ "e- 

in: 
/ = / / 



and 



^Wb - . 




°i'l. 


7£.. 

^<Cu 


Vt 


*■ £L 


^ ]} 


( 133 ) 


See list of symbols for definition of these velocity components* 
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( 139 ) 




-h 7Z t -.J +(Y~ t 


f v ^ 


*- rJ 


/ 




(ito) 
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■-) (- J(X^- +- 7 Z_ 


C,x^~ Y^ ^-J (. y^X^. ^ +'X 


Ok± Jit Zj* J 
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(i 4 i) 
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^*4 


“ ^f. / *'«Yt. 3 


r < - 

^ JCXrJ ( - ~' 7 Z OJl ) 


0 <<t- 


cz-'S 




/* 5 ~ 




-7 


(ite) 


U 9 



where - '***% XvJs ^ 


v , 

* vo. - » 


JSy^ 


XL- 


^ = i ^ *<* , £4 ” ^ ' 11 J 


% 



The indices ^ and z refer to the chordwise and spanwise source line locations , 
respectively. 

The "bracketed terms of equations l4l and 142 are seen to be identical 
to the third and fourth bracketed terms of equation l4 in Appendix C , 
resnectively. Because of this, the velocity due to thickness and the induced 
velocity due to the bound vortex segments can be computed simultaneously 
saving a considerable amount of computer time. 

It is noted that equations 139 and 136 are the same as equations 13? 
and 138 , respectively, except for 'SCtyc ^)^/ being substituted in 
for ( d %k) ) iL - Also, the last terms jn equations l4l and 142 are 
omitted when substituted into equations 137 and 138 . 


^ANPOP ON WING 


The velocity induced tangent to the wing surface by the fanpod quadri- 
lateral vortices is given by; 



where 


± (ihh) 

= />"» j f 1 ^) 


1 ^ L 

i>~ j 


(146) 

(14?) 

(143) 
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[>J 


( 1 ^ 9 ) 


since c*,.^ is small; 


{^ F } " f Vf'K^w^xig? ]{-§} (150) 


Also, 


ft?}- 




■g 




fe] 


( 151 ) 


■where 




J 


and 




‘? w \ 
Therefore : 


/ 


V; j " [ ± ]{ts} 


ir«. 


( 152 ) 

( 153 ) 


( 15*0 


( 155 ) 


The influence matrices s^h /2*r*y F are computed with equations 

derived in Appendix D. 

FUSELAGE OH WING 

The velocity induced tangent to the wing surface hy the fuselage 
quadrilateral vortices is given by: 

[%} + F T ^ ^ * 7 Cl56) 
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therefore., since 7^ = - O ‘ 


7, 




Jfej 


(157) 


Also; 


v. 


TwB ( - 


\£= 


•Tl I*J 

'* *yO *V 


7^r ~KrJ- r 


M 


(158) 


since 77 

A/— 


VfyJ 




- £> 


(%} ’ 


/ 


I C . -Lsi. \f ’j^g -f- 7^ /IiaJS 

J_Tj[i-f-Ci+t:/5^0C cijf^ - 5^9 J *nJL 


l € 


(159) 


The influence matrices AoOb and /?u 3 are comnuted with equations derived 

y - Ts 


in Appendix D- 


*3 


PYLON ON \-JlW: 

The velocity induced tangent to the vine surface hy the pylon vortex 
and source lattices is given by; 



( 160 ) 
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whore? 


[ /W r„ ] = [ Ss '^s][ 'H-l - 

r^j ■ + nik] 

[^J ' [ 5--V ] 

[*^] ■- ] - r-i «■* j 

and 

] - [ w -I * h . ] * ] 

Therefore : 



(16X) 

( 162 ) 

( 163 ) 

( 164 ) 

( 163 ) 

(166) 


and 



(167) 
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since 




o 


v. 


] 



( 168 ) 

(169) 


( 170 ) 


tnen 


= . V{i s] [{ ['*"*<][ ■*-: 


*] 


r 


( 171 ) 




-[w s |s-^ 1M) 


The influence matrices and are computed from 

equations defined in Appendices C and D. The matrices , SL/7^ j 

and are computed from equations defined in Appendix B. 

NACELLE OR WING 

The velocity induced tangent to the wing surface by the nacelle source 
frustums is given by; 







+ Lt, 

If 




(- 72 ) 
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since 


~ [ / ^xj J 


[ ' s "'5w = - P" y ~ Am &> 1 
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Also; 

{¥} ’ & 
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(17*0 

(173) 

(176) 

(177) 
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‘S/*^JGrJ r J 
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/tW/v' 


>w 



(178) 


The influence matrices and 

defined in Appendix E. 



are calculated with equations 


TTPSFSTREAK ON WING 


Expressions for the velocity tangent to the wing surface due to the 
freestream were given in reference 




r a iK 

“ oTA'j W J 


[/ /- fa)C ^0 

+■ s* 7^. 

— O'X^J ^ J 


(17?) 
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V r 
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(180} 
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W7 NG OK PYLON 


The velocity induced tangent to the nylon surface by the wing source 
and vortex lattices is given by the following expressions. 

{%} s + 


where 
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( 181 ) 


( 182 ) 

(183) 

(18*0 

(185) 

(186) 

(187) 

(188) 

(189) 
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Therefore ; 


fet 




+ IH- ]{=}} 

(191) 

Also: 

fej- 

^ jfe] 


since 

+ i‘S> s K, + 

(192) 

H 

OH" H 

(19!) 


(19*0 
(195) 

Therefore ; 

{-§rj ' [ t OfH Ife] 

jfefj d#) 

The influence matrices , /^johJy^ , and ^w ?>J are computed 'by- 

equations defined in Appendices C and D. The matrices , 

and are computed by equations in Appendix B. ft ‘ v ^ 
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FANFOD OK PYLON 


The flow induced by the fanpod quadrilateral vortices tangent to the 
pylon surface is given by the following expressions. 


{£=} = 

[V 1 '-s. / H. t ^ ^ r7 . 


(197) 

where 




f /5Wjr r 

] * [ ] 


(198) 

f/H= 

} +_ 


(199) 




Therefore ; 


[y - [' * saw* Jfel 


also; 


~ i T ?< r 

wij /_* 


r 7= + ^ /?7 ’ /r » 




(200) 


y = -1^ J(s 


Hie influence matrices i s%3f and -''*5=^ are computed from 

equations defined in Appendix D. ^ ^ P ~ roni 
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FUSELAGE ON PYLON 


The flow induced "by the fuselage quadrilateral vortices tangent to 
the pylon surface is given by equations analogous to those for the fanpod 
on the pylon. 

Therefore; 



and 






(20U) 


( 205 ) 
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The influence matrices 




r P3, 


equations defined in Appendix D. 


*3 




, and are computed from 
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PYLON ON PYLON 


The flow induced by the port and starboard pylons vortex and source 
lattices onto the starboard pylon is given by the following equations. 
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where tan?^ _ and tan are defined in Appendix I 

The slopes ; , and are determined numerically from the 

spanwise variation of the coefficients 4^ , 4^ , and ^ obtained in 

the simultaneous solution of the influence equations. 
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The expressions for ^ ^ , and ^ are the same as for 

the wing except that the subscrip?. ?=> ^""substituted f or ^ . The equivalent 
wing expressions are given in equations 135 , 136, 137, and 138, 

The matrices and are due to the port pylon and are 

computed from equations defined in Appendices C and B, respectively. 


NACELLE ON PYLON 


The flow induced on the pylon surface by the nacelles is given by the 
following expressions. 




(217) 


where 






-to 


r / H'] 


N 




[ 




Y 


+ 


(218) 

Ib'+J (219) 


and 

r 




r. 




"\ 




( 220 ) 


where /9-^u and ^ 4. are defined by equations (I 9 ) and (50), respectivelv. 

r 

therefore; 
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then 
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The matrices /fatL and are computed from equations defined in 

Appendix E. ^ ^ 
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FREESTREAM OH PYLON 


The expressions for the velocity tangent to the pylon surface due to the 
freestream are the same as for the velocity tangent to the wing surface due 
to the freestream except for the substitution of the subscript p for w . 


+ ± ) * J yz 

“sr _ p ***** y^X' + 
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( 224 ) 
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WING ON FANPOD 

The velocity induced by the wing source and vortex lattices tangent to 
the fanpod surface is given by the following expressions. 
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Since the 'ding and fanpod coordinate frames are oriented the same 
direction; 



(228) 


and 



The influence matrices , and are computed 

•using equations defined in Appendices C and D. Also; the influence matrices 
3 > an ^- are computed using equations defined in 

Appendix 33. 

The unit vectors and tangent to the actual fanpod surface 

are determined from the unit vectors and } tangent to the fanpod 

surface obtained by the Gothert transformation, in the following manner. 
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™ere -— — ' and - . - - are obtained as follows; 
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Then from equations 230 , 231 , 232 , and 236; 
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and after substituting equations 233, 23b and 235 into equation 238. 
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FAHPOD OK FANPOD 


The velocity induced by the fanpod quadrilateral vortices t ang ent to 
the surface of the starboard fanpod is given by the following expressions. 
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(aid) 


where and are the values of the surface vorticity along and 

perpendicular to the meridian lines, respectively. Hie vorticity is determined 
by dividing the average of the two adjacent vortex strengths by the average 
perpendicular distance between the two adjacent vortices. 

The influence matrices /f/syr > > and are computed with 

equations defined in Appendix D. The components of the unit tangent 
vectors and % r are computed from the vortex grid coordinates. 


PYLON ON FANPOD 


The velocity induced by the pylons tangent to the fanpod surface is 
given by; 
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and 
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The influence matrices , and >3=7^_,are computed with 

equations from Appendices C and D. Also; the influence matrices ? 

S-P 3 , and Si— are computed with equations from Appendix B. 

% ' “^p 

NACELLE ON FANPOD 

The velocity induced by the nacelle tangent to the fanpod surface is 
given by; 


~ ' if- + i rp "v F ^ ^ 


(250) 
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The influence matrices and /%=^y are computed with equations 

from Appendix E. ^ " y ' J 
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PRESS TP ELAM C” FARC'D 


The components of one 
the fanpod are given by; 


freestream velocity tangent to the surface cf 



-V- 




( 255 ) 



(2 56) 


The veloci 
;o the fuselane 


y induced by the wing source and vortex lattices 
surface is given cy the following expressions 
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The influence matrices , and are computed using 

equations cefinea in Appendices 0 ana 0. Also the influence matrices <5^ 

■ and ^ w #«, are coni P uted using equations defined in Appendix B. X " ’ 


The unit vecotrs ( , ~ 1 T3 r ^ r , 

l-p^ ) are computed in the same manner 
"*T3 


T S%) and ( ' 7 ^t Aw ; 
as those for the ?anPod« 
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FAI'fPOP ON FUSELAGE 


The velocity induced by the fanpcd quadrilateral vortices tangent to 
the surface of the fuselage is given by the following equations. 



and 

[%} 



(261) 


(262) 


where the influence matrices - y , and are computed 

with equations derived in Appendix D. r k '" A 


FUSELAGE ON FUSS FA OE 

The velocity induced by the fuselage quadrilateral vortices tangent 
to the surface of the fuselage is given by the following expressions 
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Where and ^ 5 y. are the values of the surface vorticity along and 

perpendicular to the meridian lines, respectively. The vorticity is deter- 
mined by dividing the average of the two adjacent vortex strengths by the 
average perpendicular distance between the two adjacent vortices. 


The influence matrices ^•33 
equations defined in Appendix D. 
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y , and 


are computed with 
*5 
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PYLON ON FUSELAGE 


The velocity induced by the pylons tangent to the fuselage surface 
is given by; 


[%} ' [' 


7 i T ^ x 
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( 2 66 ) 


( 267 ) 

(268) 

(269) 

( 270 ) 

( 271 ) 

(272) 


The influence matrices , and are computed with 

equations from Appendices C and D. Also, the influence matrices , 

, and 'S^-p> are computed with equations from Appendix B. ** 

Yt> * 7 > 
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NACELLE ON FUSELAGE 


The velocity induced by the nacelle tangent to the fuselage surface 
is given by; 

+ r/-r B + 


( 273 ) 


and 




+ i + i X/H J 
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where 
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The influence matrices ^3^ and are computed with equations 

from Aopendix E. ^ jUj 


FREESTRFAM 0W FUSELAGE 

The components of the freestream velocity tangent to the surface of 
the fanpod are given by: 
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and 
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COMBINED WING VELOCITY EQUATIONS 

When the velocity induced by all the configuration components is summed; 
the following expressions are obtained. 

In the cnordvise direction, 




( 280 ) 


and in the spanvise direction 




J 


(281) 



are defined by; 


Rr} = [l.K.w^xHr^ 




and 






d £7 J 


(282) 
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Then, the velocity components in the cbordwise and spanvise directions, 
corrected for compressibility by the method of Labru,iere, are given by; 
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COM BirrSD PYLON VELOCITY EQUATIONS 

The sum of the velocity, induced by all of the configuration components, 
or. the nylon is given by: 


(%] *{^3 ovi + [% , j 


(290) 


in the chorawise direction, and 
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( 293 ) 


then, the velocity components in the chordwise and spanwise directions, 
corrected for compressibility by the method of Labrujere, are given by; 
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AV„ r (295) 

* O* ¥■ i%, j 

where ^ is defined the same as on the wing except that fa is used 
instead of • 

COMBINED FAKPOD VELOCITY EQUATION 

The sum of t^e velocity induced on the fanpod surface by aJ 1 of the 
components is given by; 

(%} + fe] f f%] + f W + (296) 


alo ig the meridian lines and 

«- [%] + ^ 

perpendicular to the meridian lines. 

COMBINED FUSELAGE VELOCITY EQUATIONS 

The sum of the velocity induced on the fuselage surface "by all of the 
components is given by; 
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along the meridian lines and 



=[%] + (%1 + [%} 


J 


(299) 


perpendicular to the meridian lines. 


The compressibility correction for the fuselage and fanned is that 
due to Gothert and is already included in the above equations. Due to the 
assumption of axisymmetric flow for the nacelle, no velocities or pressures 
are computed on the nacelle. 


The surface pressure coefficient is then computed at each of the 
control point locations using the following expression. 


Ycv-o 
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(300) 


where ^ 


is the ratio of specific heats. 
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SECTION AND TOTAL LOADS AND MOMENTS 


The section loads on the fuselage and fanpods are computed by summing 
the products of the surface pressures and directed incremental areas at each 
longitudinal station. 


The section drag is given by; 

(%£) = - ^ 


( 301 ) 


and the section lift by: 
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where c is indexed over all subareas at the longitudinal station k, and 


a = M'SJ 

te. i 

is the fuselage or fanpod reference area. 


( 303 ) 


^AVS- // ^- 


The fuselage or fanpod total lift, drag, and pitching moment are 
computed as follows. 
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r r 1 he section loads and moments on the wing and pylons are computed by 
the following equations at a series of span stations* 

The section lift is given b:/; 


Cn G c /9JS f 77 ' 45 , , 


"A/G 

the section drag by; 
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and the section moment by; 
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where 
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::CKEt y is equal to 90 degrees on the wing. 

The section center of pressure is given by; 
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Q jJL 
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( 312 ) 


The wing and pylon total loads and moments are computed by numerical 
integration of the section loads. 

In addition to the drag computed from equations 301 and 309 , the 
program also computes the skin friction and vortex drags for the complete 
configuration. The equations for these calculations are given in 
Appendices J and K. 
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EXPERIMENTAL COMPARISONS 


Two configurations were run to demonstrate successful program operation: 
Case I - The wing-body configuration shown in figure 7 was run at 0.6 Mach 
number and zero and four degrees angle of attack. Case I comparisons are 
shown on figures 8 through 27* Case II - The v/STOL lift fan transport 
model (less tail surfaces) shown in figure 1 on page 13 was run at 0.6 Mach 
number and at zero and three degrees angle of attack with and without the 
fanpoa. Case II comparisons are shown on figures 28 through 57. 


WING-BODY CONFIGURATION 

The wing-body configuration run as Case I is defined in RM L51F07. 
Comparison of results from the program with data are shown in figures 8 
through 27. Figures 8 through 12 show comparisons of pressure coefficients 
at zero degrees angle of attack at the 20, 40, 60, 80, and 95 percent semi- 
span wing stations. 

Figures 13 through 16 show both the body alone and wing-body pressure 
coefficients at the 45, 75> 105, and 135 body roll stations at zero angle of 
attack. These figures show the increment in pressure induced on the fuselage 
by the wing thickness. 

Figures 17 through 21 compare results of the program with data at the 
20, 4o, 60, 80, and 95 percent wing semi-span stations for four degrees 
angle of attack. Figures 22 through 25 show both the body alone and wing- 
body pressures at the 45, 75, 105, and .135 body. roll stations at four degrees 
angle of attack. These figures show the pressure induced on the body due to 
thickness and lift of the wing. 

Figure 26 shows the comparison of theoretical to experimental span load. 
Figure 27 shows the comparison of theoretical to experimental longitudinal 
loading on the fuselage with the wing at four degrees angle of attack. 

All of the results for Case I are considered to be in excellent agree- 
ment with the test data. 
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MAXIMUM DIAMETER 


Figure 7 k Wind Tunnel Model 
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FIGURE 9. CHGRDWTSE PRESSURE DISTRIBUTION AT THE WING 40 PERCENT 
SEMI-SPAN STATION- * ,6, OC« 0 DEGREES 
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FIGURE 11. CHORDWISE PRESSURE DISTRIBUTION AT THE WING 80 PERCENT 
SEMI-SPAN STATION. « .6, oQ = 0 DEGREES 
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FIGURE 12. CHORDWISE PRESSURE DISTRIBUTION AT THE WING 95 PERCENT 
SEMI-SPAN STATION. .6, oO 0 DEGREES 
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FIGURE 19. CHORDWISE 
SEMI -SPAN 
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PRESSURE DISTRIBUTION AT THE WING 60 PERCENT 
STATION. = .6, cC- 4 DEGREES 



















TEST 

DATA 


O 

□ 


PRESSURE 

STATION. 


DISTRTBI.m ON AT THE VTN'O 95 PERCENT 
« .6, oC* ^ DEGREES 























AT THE FUSELAGE 75 
, OC« u DEGREES 






JIJTTOH 
* - 6 > 


AT THE F> 'SETAG 
OC= h DEGREE:3 











PRESSURE DISTRIBUTION AT THE FUSELAGE 135 


L STATION. M » .6, CO 4 DEGREES 


















FIGURE 27. 
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V't-TOL LIFT FAN TRAIfgPQRT 

The wing-body and wing -body- fan pod configurations run no Case IT are 
-"fined In NA-72-17-2.* Comparisons of pressure coefficients at zero and 
■fiee degrees angle of attack were made for the wing-body, the wing-body- 
fanpod, and the difference between the wing-body- fanpou and the wing-body. 

Figures 23 through 32 show comparisons between data and the program at 
the 20, 33, 60, 75, and 95 percent semi-span wing stations, at zero degrees 
angle of attack, for the wing-body configuration. Figures 33 through 37 
snow analogous comparisons for the three degree angle of attack case. 

Figures 38 through 42 and figures 43 through 47 show comparisons between 
data and the program at the 20, 33, 60, 75, and 95 percent semi-span wing 
stations for the wing-fuselage -fanpod configuration at zero and three degrees 
angle of attack, respectively. The 33 and 60 percent semi-span wing stations 
are just inboard and outboard of the fanpod, respectively. 

The comparisons for the upper surface are quite improved over that for 
the previous program. However, the lower surface comparisons are still not 
in very good agreement. This difference could be due to either thickness 
effects of a high wing configuration, which can only be properly accounted 
for by a surface singularity approach, or due to a wing-fuselage juncture 
interference problem in the analysis. It is believed to be the former of 
these two situations, because the problem does not occur in Case I which is 
a mid-wing configuration. 

The increments in the wing pressures due to the fanpod are shown in 
figures 48 through 52 for the zero degree angle of attack case and in figures 
53 through 57 for the three degree angle of attack case. These comparisons 
are not as good as expected. The three degree angle of attack case is 
considerably poorer than the zero degree angle of attack cose. Also, the 
comparisons are worse right at the junctures. Both of these situations could 
be indicating boundary layer effects, since the boundary layer would be 
thicker at the higher angle of attack and right in the juncture. 
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FIGURE 32. CHORDWISE PRESSURE DISTRIBUTION AT THE WING 95 PERCENT 
SEMI-SPAN STATION. .6, oQ* 0 DEGREES 
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FIGURE 33. CHGFDWISE PRESSURE DISTRIBUTION AT THE WING 20 PERCENT 
SEMI-SPAN STATION. M ^ * .6, oC » 3 DEGREES 
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FIGURE 3 It-. CHORRFISE PRESSURE RTSTRIBUTION AT THE WING 33 PERCENT 
SEMT-SPAN STATION - . M ^ * A, oC » 3 REGRESS 













Ill 



FIGURE 37 . CHORDWISE PRESSURE DISTRIBUTION AT THE WING 95 PERCENT' 
SEMI -SPAN STATION. M w » .6, oC * 3 DEGREES 
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FIGURE U 3 . CHORDWISE PRESSURE DISTRIBUTION AT THE WING 20 PERCENT 
SEMI-SPAN STATION. M ^ ,6, o<T * 3 DEGREES 
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FIGURE U8. PRESSURE INDUCED ON THE VI NG BY THY FANPOD AT THE WING 
20 PERCENT SEMI-SPAN STATION. M ro = .6, oO 0 DEGREES 
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FIGURE 51. PRESSURE INDUCED ON THE WING BY THE FAN POD AT THE WING 
75 PERCENT SEMI-SPAN STATION. .6 qO* 0 DEGREES 











FIGURE 53. PRESSURE INDUCED ON THE WING BY THE FANPOD AT THE WING 
20 PERCENT SEMI -SPAN STATION. .6,0c* 3 DEGRESS 









FIGURE 5 I+. PRESSURE INDUCED ON THE WING BY THE FANPOD AT THE WING 
33 PERCENT SEMI -3 PAN STATION. K ^ *= .6, OC- 3 DEGREES 








129 








FIGURE 56. PRESSURE INDUCED ON THE WING BY THE FANPOD AT THE WING 
75 PERCENT SEMI -SPAN STATION. ■= .6, oC x 3 DEGREES 









FIGURE 57- PRESSURE INDUCE" ON THE WING BY THE FANPOD AT THE WING 
95 PERCENT SEMI-SPAN STATION. M C?Q = . 6 , 0 C~ 3 DEGREES 






APPENDIX A. 


NUMERICAL PROCEDURES 

Discussions of the prime numerical procedures used within the program 
are given in this appendix. There are essentially three such procedures ; 

(l) straight line interpolation and extrapolation, (2) controlled deviation 
interpolation, and (3) Householder’s simultaneous equation solution. 

For straight line interpolation and extrapolation about two given points 
(Xi, Yi) and (X 2 , Y 2 ); 


y-=. cl +- i ~ ; Y, 


( 1 ) 


where 

X— X, 

* & x*-x, 

The slope for this case is given by; 

57 - *'CY x -r,) 


(a) 


(3) 


where 

■*' = VOk-x,) (4) 

In the case of the controlled deviation interpolation method (CODIM) 
parabolae are used to curve fit a set of four given points (Xjj-1, YN-l), 

(X N Y W ), (X K+1 , Y K+ i), and (Xn+2, Yn+ 2) to obtain interpolated Y and ^ 
values f t r Xjj £ Xg^ +1 * Only that information, relative to this method, 
which is necessary to judiciously pick input points will be discussed here. 
A complete derivation is given in reference (20). 


One parabola Pj is fit through (Xn_i» Ypj-l), (Xn, Yra), and (Xn+ 1, YN+l) • 
The other parabola Pq is fit through (Xn,Yn) , Xjj+i, Yn+i), and (Xjj+2, Yn+ 2) • 
This curve fitting process involves the solution of two sets of simultaneous 
equations . If, 

JP as. /4 f x ' 1 Y- X -f- C f (5) 


and 

7| = x 1 /- 
Then 
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a nr? 



The interpolated values of Y and 5^. "between Xj,j and Xfj+l are defined by 
either Pi, P 2 , or a linear combination of F^ and P2. The amount of Pj and F2 
used in the linear combination is determined by comparing both of the 
parabolae to the straight line. 


S - <* Y^-j- 

r + Cl-=t; 

( 9 ) 

where 

0 c = 

XpJj- f 


(10) 


The parabola which has the least deviation from the straight line is 
given the greatest weight. The weighting factors Ex and E2 are determined 
as follows; 


- /?-*■/ 

(11) 

^ = /?-*/ 

(12) 


The weighted expression for Y in the range Xn ^ a is then; 
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The derivative ~~ in the range X$j ^ x ^ is then; 

c?Y cf/o / , d~p / -y 

" J^/td “ aJ 3 H/' d *** 


( 13 ) 


( 1 *) 


where 

/^’= 7f -f- / - ^ J ^ ^ 

•&=<*&, T 4 C ' ~ 

And then 

dX dX'^- d* ^ ^ 


( 15 ) 

(16) 


<dP 


*x zp- ? ( 17 > 
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Therefore ; 


d~D _ c/d. 

dx ~ dx ' 


<dX 


c/ai 

3jr e * 




(18) 


r, 


y_ J «t ^ i +■ O -°0 "p, 

°i s t +- C i -«0 

V- TL 


Fotz x-x^j 


X&7Z X^j < x < x^, 




X - Xfj+ r 


(19) 


and 



dr° / „ , ctr> / 

3*A> ^ Jx h 


FC7Z X^Ssj 

x^c x < X^4. , 


F07Z X=X^ / 


( 20 ) 


In the case of an end interval Xjj_x ^ < ^rO > P 1 is se-t e( l ua l j 




( 21 ) 


where 

//*%/ -/^* / / 

A = / - 

/Si,/ *-/***/ 

and 


~~ Kv~/ 


~~ XrJ-l 


( 22 ) 


(23) 


YrJ YrJ*-/ 


“ Xj^-x 


rj+ t 


(24) 


A similar procedure is followed for the other end interval X N+ x^ X< 
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Householder's method for solving simultaneous equations is used in the 
solution of the aerodynamic influence equations. The method is applicable 
to both square and rectangular influence matrices. In the case of 
rectangular matrices it Is not necessary to least square the equations first, 
since Householder's procedure least squares and triangularises simultaneously. 
Also, the influence matrix is triangularized by means of orthogonal 
transformation matrices, which preserve the conditioning of the matrix. The 
combination of these two advantages, along with a reduction in the number of 
required computer operations, greatly improves the numerical accuracy and 
stability of the solution over that of the standard Gaussian reduction method. 

A complete, but rather abstract, derivation of the method is given in 
reference (21 ) . The actual computer subroutine used in the program was 
developed by L. V. Andrew at NR. The method In the subroutine has been 
altered from the original to allow the operation on a single row of the 
matrix at a time. This reduces the required core allocation necessary to 
triangularize the matrix. 

A derivation of the method, developed by the writer, will be given here 
in order to described the basic philosophy of the method. 

If Is the rectangular influence matrix, the upper triangle is given 

by; 

1>7 (25) 

where is the combined orthogonal transformation matrix used by 

Householder to triangularize £ A J . 

The relationship between Householder’s triangularized matrix ( iz ) 
and that obtained by Gaussian elimination of the least squared influence 
matrix £ A } , is 

D 3 -] = t> x* ft* 1 ~ £>*7Z> 72X7 =7^ 7&J ( 2 6 ) 

where is the triangular matrix obtained by Gaussian Elimination of 

[a j[a 3* The matrix is the Gaussian transformation matrix used to 
triangularized ] [A !• And, X} is a diagonal matrix with the same 
diagonal as [72.}, It can be seen from equation (22) that a nonsquare matrix 
must be least first, before applying the Gaussian transformation Ct], 

Where as, the Householder transformation matrix can be applied 

directly. The least squared matrix IaYI*} is usually more ill-conditioned 
~:*an (70, and therefore, less accurate results are obtained. 

In the Householder method X3 is equal to the product of N+l Individual 
orthogonal transformation matrices, where n equals the number of unkowns. 

There are N+l transformations because the augmented influence matrix, made up 
of the influence matrix itself plus the boundary condition matrix, added on 
as the last column, has N+l columns. Each transformation results in reducing 
all elements below the diagonal to zero for one column. The columns are 
reduced from left to right. 
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The individual transformation matrices are defined "by; 


o - ct* i - -z t ; (27) 

where is a unit diagonal matrix and is a column matrix defined hy 

the unit vector^*<^“ £ C7^^. The vector 5^, is defined by the mth column 
of [a ] where the elements on rows less than m are replaced by zeros. The 
unit vector t£, is defined by a column matrix {V} w with all zeros except 
for the mth. row, which is equal to one* The constants and are 

defined as, 

(28) 

“ l/ z ' TC, (29) 


It can be shown that is reduced to /f<Z~*3 / is premultiplied 

Also, that the first m-1 rows of 

z>~. x 1 7 • • • z>, 7z>tj (30) 


remain unchanged by the mth transformation. The result after m 
transformations is then zeros below the diagonal for the first m columns 
and /£&, 7/, /(a.\ ^ /, on the diagonal. The 

elements above the diagonal have been defined by the m preceding 
transformations and will remain unchanged for the N+l-m remaining 
transformations . 

(<3x7 - 7 V-J (31) 

47, “ <"47, ^)/ y ^ (32) 


^ = (33) 


and 


/i, " aae (3^) 
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remain to be proved. It is helpful in tne derivation of equation ( 31) if 
tne vector identity 

?C*trr> =4^7 (35) 

is observed from the following vector diagram. 





Then from equation (35); 


+ Z (36) 

y 2 4^ (37) 


Tneref ore ; 




(38) 


where and are dyadics. The unit vector C^r, is in the 

direction of <2 m . 

Equation ( 38 ) can then be written in matrix notation as follows; 


O-ij- * [*l£* J£,)fcL - 


(39) 


Widen is equal to equation ( 3-1) • In matrix or tensor notation it becomes 
evident that the dimensions of fa 7 /V 7 and /«? are not limited to three. 

( 1 + 0 ) 

and 
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(4-2) 


Then if equation ( 31 ) is premultiplied by J T 


And substituting =w? and 

^ - JsvjL = *■£,, 
or 

" yf C°Cm ~ J 

In vector notation equation 

" J °^*TT1 C^st*^ ” ‘&srr7 


into equation (42). 


( 4*0 is seen to be equal to; 


(43) 

m 


(45) 


Also,, if equation (44) is substituted back into equation (31) 

~ "ftzSjfir, }) [«—] - [' 2 C~J (46) 


Therefore; 


"J °VaT-7 

J? <<** (<*„, - fa~,} 

or in vector notation 

_i» ~ ^rrT) ~^>-7 

(°trrt ,S ?*" ^ 



(47) 


(46) 
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APPENDIX B 


SKEWED SOURCE INFLUENCE EQUATIONS 


finite source lint: 


The potential for a finite source line can be obtained from the 
integral of the potential produced by a distribution of three dimensional 
sources placed along an axis. _ An expression for the potential in an incom- 
pressible flow field will be derived first, the partial derivatives taken 
to obtain the perturbation velocity, and then the velocity expressions 
transformed by the C-othert similarity rule to obtain subcritical compressible 
velocity expressions. The potential for a three dimensional source Is given 
by: 


'O'* c %- < 


n 




( 1 ) 


where f ^ Y g } £* ^ and ( j Y s j % ) are the locations of the point being 
influenced and the source, respectively. The prime indicates the local 
coordinate frame shown in figure B-l. If a distribution of constant 
strength sources is placed along the x' axis, the potential of these sources 
represent that of a finite source line and is given by the following integral 
expression. 


Xs 


cUj 


[(Xg-Kf + %■*+ 


(a) 


where (TCX S ) is the source strength per unit length along the x' axis from 
X " -S to A''- S . 

Let and h ~ Y g % , then for b ^ O 




(TCX') 

47r 




( 3 ) 
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The perturbation velocity parallel to and those ‘^"'and perpen- 

dicular to the source line are obtained by talcing the partial derivative of 
equation ( 3 ) with respect to X^ , y' } and , respectively. 

f / gT * 


u 


47r l 


+ r*+ f 7 * 


CM 


^ roQf VRK-*?* Z?+ €J' A 
47r l ( *>-*) + lcx $ - sU f V 


_ %Z£tK+*f+ K'^ % 7 * 1 


( 5 ) 


• rcxl) f %tJ &5rfl * % * jl Z_ 
s ^ C <*i- s >+t<*f s 5+ ? v ?7*‘ 


eVfcy^V a;7* ? 

<^ s ; +£cj(**f+ $'V «;7 4 J 


(6) 


The flow produced by a finite source line is axisymmetriCj therefore, 
equations (^), (5), and (6) can be reduced to just two velocity components. 
Tlie component u* parallel to and that perpendicular to the source 
line are then defined by: 

f rCXs) / _ ^ 

* = ~4¥j? ( ces % ~ cos ?;) (7) 


* 


I 


rcx s } 

4-TTh 





( 8 ) 


1 U 0 


■where <9 and B are defined in figure 1. 

e 1 f- 



Figure 1. Finite Source Line 


Equation (8) is seen to "be the same as that derived from the Biot-Savart 
relation in reference 8 for a finite vortex segment. The influence 
functions for the velocity ’V'* , due to a finite source line, and that for 
the velocity perpendicular to the plane containing and , due to a 
finite vortex segment, are therefore identical. 


CRAKKED SOURCE LIKE 


In general the source lines vill he cranked to follow the constant 
percent chord lines of the wing, as shown in figure 2. since this is a 
planar analysis, all of the source lines will lie on the chordal plane. 

The strength of each source line segment is equal to the integral of the 
distributed source density over the panel it represents. Each source line 
segment is divided into two spanvise sections; the inboard or port half 
denoted by c and the outboard or starboard half denoted by <3 , If there 

is a crank in the source line it occ\irs at the midpoint of the segment. The 
panels or source line segments are equally spaced in the spanwise direction. 
In the chordvise direction the lines are placed at equal increments of O , 
where © - cos' (/- zx/c) 

The expressions for the velocity induced by a single cranked source 
segment due to the inboard half ( u* “ir * ) and the outboard half ( «-t_, VS * ) 
are obtained directly from equations (7) and (8) and figure 3 . 
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¥■ 





Y 




V * 



Figure 3* Cranked. Source Line Segment 




/ 

J(x r x s mY^ Yj+c%- 


/ 

jR*;-# + yfwjj'+Z y^k+zF+z%-*x 



(9) 


1^3 



/ 


*4* 


<ra, %x 

4rr 


{jLf-h- v - 1 - <A* ** ' *>' 


/ 



Equations ( 9 ) and (10) wore obtained by noting that Equations (^) 
and (7) reduce to: 




l = “ 


<rc.Q 

+7r 



Also; 



4rr 




JLcXfO -cy Ctf 


{ 


Yz~Y p +<x,-x s )r**>& + %Q* r**L>*#J 

***>$/'+ £Y S - Y a + 3 J'+tY'ZJ* 


K~ Y* * < x,~ 

J~Of * a f+ c %~kx *- < % - 



* rc^Ysi 

*a — • 


4Tr sjB-XfXs) ***?%) 


{ 


%- y, + e*, - xj f c 




_ Y,~ Y, *‘Cx r -K)r*'Vfj- vC/s sw^it) 7 

'jRXfW- %.r#«j$''+ZY-y r yj'+r%-tl z J 


The following values of K and were substituted into equation (8) 
to obtain equations (12) and ( 13 ). 


A * sv * <% 


(no 


K = Jr<Xf - cr s - + C? - ?/ 


(15) 


The velocity components in the prime coordinate frame shown in figure 
are obtained from the following relations. 


/ f 

4f 

(16) 

vr'=r irc&'S y 

(IT) 


( 18 ) 

Therefore ; 


< = 

(19) 

- V(T<1 % -v; 4 - C r f ~ xj_7 

(20) 

^ " < £- if *~C\-?JO+ r*Aj%) 

(21) 

, ~*rt C X - - c* s 

•» VE111. iijj _. .. i i ■ 



(22) 

•zrV ?r - a;-/ ' *• 

‘ v£Oy-^ - c y~#c*«>?J'+ c %- 

(23) 


, irk $ - ?)JT+ 

w °~k£Cy s - XJ - ly- v ( ? ? - zjc * + 


( 24 ) 



Then using the following transformation equations the velocity component 
in the freestream direction ^ , that in the spanvise direction ’ir , 

and that perpendicular to the chordal plane w are obtained for a single 
cranked source line segment 


CL ~ 




cejZsPsJ -- -Tf ' 

^ V * + ~ 


(25) 


•2r - 


i/. + -v? r s? A) & 

y / ^ tv?*?# 


/~77- zrs?s^£ 


(26) 


w t= w'- + uvj 


(2T) 


SOURCE STRENGTH 

As mentioned above the source strength in a planar analysis is only a 
function of the local symmetrical airfoil surface gradient. Th: s is due to 
the fact that sources cannot produce a velocity component perpendicular to 
the chordal plane at locations oth '.r than their own. Also, in order to 
satisfy continuity the source strength must be equal to the mass flux emitted. 
Therefore, the local chord wise source density is; 


W*jK) = ^ r 5 j 


( 28 ) 


The discrete source is assumed capable of representing a linearly 
varying source distribution in the same way that a discrete vortex represents 
a linearly varying vorticity distribution. Therefore, let 
be defined by the first two terms of a Taylor series expansion about the 
location of the discrete source. 




^ X 


(29) 


"Where the subscript s on <r indicates lines of constant source density 
parallel to the discrete source line. 


Also, by equation (28) 

= ^C^Y S ) ^ LX~**> 


(30) 
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The discrete source strength Yc) is obtained "by integrating 

equation (29) over the panel in the chordvise direction. 


<n^Y s ; = 


= f (. X J 


$S “*• \ V. ^ J yjcos cf> d a 

o 


= j"r<s<^s; + cx-^ ] C c~ ? jk 


(31) 


where Jl - 5 ' ^ &s is the panel chordvise length and © 5 ~ cos 0" 

The qos j25 terra appears due to the fact that the integration is being taken 
in the chordvise direction rather than perpendicular to the discrete source 
line. 


- £<£<>,, Y S )JP 


z — 1 M 


(32) 


since and Yj = Z j 


Y / 7^ ^ 


(33) 


If the source strength is divided by the free stream velocity and the 
boundary condition for no flow through the airfoil surface ~ ^~-C^ Sj Y s ) 

is substituted into equation (33) the following expression for source 
strength per unit free stream velocity is obtained. 


7^V,) =■ 




o'? 


s^c/7+ 




( 34 ) 


where equals the number of source lines per chord. 
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COMPRESSIBLE VELOCITY RATIO EQUATIONS 


The final expressions for the compressible perturbation velocity ratios 
are obtained by transforming the coordinates and perturbation velocities by 
the Gothcrt similarity rule and by substituting equations (9), (l0)> (~ 2 ); 
(13) ; and (3*0 into equations (l<?) through (27). Also, the contributions 
to the velocity ratios at an influence point ( ; ^ ) by all the source 

line segments at ( Y s J ) u are algebraically summed for = ^ ^ /wj 

and 4 rs/jg, • . 


M M 


^ ^ f, , t- 

W'skrfeZ* & *&<**?)] - <• >[ 


js g, V yj - y, V 


g£ji ^ Sk jj fcj. ^ £ C K %) 1 

J 


(35; 


Al /v< 






£>■<- A g, i 4 U o *• £ I, £) 


- ^ K * Y ’^ V & g/ jo ^ „ Kj y 

?< x ’o Yj 


(36) 


148 


and 






^,<X U Y,J 


?<X XJ 


(37) 


vhere: V Vf % , 9^%/^ > % = %/ % » 2,,= Vfe* > 

K.“ K/% , ?=f°*-v4 > X*4 OM fi ’ ? = p * anfl 

? +(Z2££ff'* 


Also; define A “ J( g - X, t . , Y A ~ % ~ ^ > and - 


5 = 


dlien for ? <3 

fl 


(*ts *2 j *0 


M 


x, + V a J z + ( K + 'J z * ^ 


y w~^ z*+ f 


(38) 


li+9 



&U J. s ^2 J = 


/(X; ~ CZ-O' f 


■J // ^ r‘ * s’ 


(39) 


<%•<*. s; = a >*l 

6 ‘ * r> s .t _ -t 


/(*<.•“* / ? ^JJcJ in^JV (y f * ? z 




< V * ^/c*. - j?^/ * flSO^v £ v e 


v, i- 


C^ t - UOJft*- 9J + * 


(*» 


%j %j **. j K.) =■ (f ^ ^ X £ »*• <C^ ~ K4j^ 

* ^ 1 F/ 7 r _. . l -i 1 A- ' ".t- • — "71 ) Sf 


Z^" ?* 4 / * *\J/(X< -Z*f + CK~') X + is* 




- YJ*j£<-*r ?£j+ eV//*>>v a - 


Vi ** <*i 73jf, 


< £ - n v ?J+ C 


(M) 
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S %. (K + KXu) + 5 ^ 


[c*<- yj<J+ *■<% + *>'+ * 


g 3 


<-X‘-ZZj{cA~K* t /+ *>IU*- + F+ ?- 


?c x «* *.-%< ; ?g. 


<^K^/Av gW 


(42) 



Cx. ~ 9 r. ? + ?V7.AJ?. 





(< / ? 3 


?c>v ■sX/yWTWr 


C K ^ •*< ^°it. ) 


<A - g ^ fwn^r^ 


(43) 
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As discussed in reference 2 the form of the influence equations (4o) 
through (43) are such that cannot he set equal to zero to obtain the flow 
on the surface of the wing without losing a substantial portion of the 
solution. This problem is due to the singular nature of these equations. 
Equations (42) and (43) will not be discussed further here since 

The first and third terms of equations (40) and (4l) cause no special problem 
as ? — *- o . However, the second terms in these equations do contribute 
nonzero finite values to , a) and Yg j eO at the tips 

and at span stations where kinks or cranks occur in the source lines. It is 
suggested in reference 2 that the influence equations he integrated in the 
free stream direction before letting t*, — **■ o . This integration was discussed 
in detail in reference 2 for the semi -infinite sheared source lines with 
~ The_same procedure used there to account for kink, crank, and 

tip effects when ^ * a will also be used here. 

If equations (38) through (4l) are specialized to the case of = © 
and the second terms in equations (40) and (4l) are replaced by expressions 
derived in references 2 and 6 , the following expressions are obtained 
x or S u . C ,x g , y 3 t Oj x lcJ 

%> Oj K; K) . 


) , and 


%j 0 j K.J ~ 


fc*i + V c k* tf /*7 + ?: 


(44) 


/ 0 j X’M Kj.) — 


c K -if /# + K 


m 


i ^ j 0 J ~ 


/ r 




- K%J l JCK-f Z,S ^ CZ + O' JJT+ K 


zn VC-cr) [,oc,s c * /t-jxj] 


(46) 
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^ J = . * f <K * «.%.■> - g.* _ C* »*■ *2,J 7 

C x " l /^v r j 

/ g. y. (■ tt s C- ./zz /ev.Z/) i *>%/ %'. - z:_/ 

2 ^^ -C^^-j [■ “ /? %?j] 


0 * 7 ) 


The third terms in equations (46) and (47) are only evaluated when 
/‘S%?J ^ .STS' 7 C L , when < sj , where -i is 

the number of the panel aft of the leading edge on which the velocity is 
being computed, and when is at a span station where a break in a source 
line occurs. When these conditions are not all met the third terms are set 
equal to zero. 

The expressions for the velocity components given in equations (35), 136), 
and (37) can he represented by the product of source influence matrices and 
the source strengths. I 
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Where the elements of 5* , ■sy , and are computed from 




^%fr 


4 = 1 „ «• - e. 


]~f- 


*a,- 7 “ £■% 
??■ 




( 51 ) 




(r- 


/S i 4 


7 - [— 


? ?. 


*J] 


( 52 ) 


and 

S ~ = ~ ^ 53 ) 
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APPENDIX G 


skewed vortex influence equations 

The velocity induced "by a skewed vortex in the plane of the wing can 
be derived from the Biot-Savart law. The expressions for the velocity 
relative to the wing coordinate frame will be derived for the incompressible 
case first and then modified by means of the Goethert transformation to 
account for compressibility * For the general skewed planar vortex shown 
in figure C-ithe total velocity induced at point ^ is the vector sum 

of the velocity induced by the four vortex segments *?<*> . 


A 







Figure C-l. Skewed. Planar Vortex 
The Biot-Savart law is given by the following expression: 





& of 'e 


CD 
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where /fij is the perpendicular distance from the vortex segment to the point 
Q, and E is the angle between the vortex segment and the line from the 
point Q to the increment of vortex segment being integrated . The integral 
is evaluated around the complete closed path co A 1 V A . The direction 
of the induced velocity is given by the following unit vector. 



The contribution from segment e^A’ is derived as follows. 

Jw = [c°sc“> - C0S-C'S&- « (3) 


where 

X — — ' Xy> 

**= V ^ 

since 
cos Co > = / 




COSVBo-ev*'#) - _ ccstoos?‘c?> - _ X* +&*&***& 

fi* 4 * Sb***#? x f +■ v*y* 

where is the sweep of the port side of the skewed vortex, 
and 


✓v 


then : 




? _ r ** •&» i 


xL 




( 5 ) 
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The contribution from segment A’V is determined similarly. 



?* (rsz-J*!*- > / >]*■ 



<? G>-s(/a& - £p iv? ^ 






( 6 ) 


where 




y* i' -a y-^/?so 4 - 
^9/jV’ J V <ry"V y^j 2 ^ 


£ CPS>O 3<0 - £?\//?'J — 


y* / x*r/^y=*- 
2?^^. +*J'*fi<**+ y* V 


and 


V; " 


2?^* f'j] 


7 T ;?*C - ? -(Tx^y^^J 


i-)*Sj 
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therefore; 



<? 4 X **>**** i-J 




g!± ***»^ 4 Iffy 

^'Vv j-yt/Z * 


&J ~ (X*- 


(T) 


The expression for the contribution from segment VA is analogous to 
that of A’V. The angle fa is the sweep of the starboard portion of 
the skewed vortex. 




y* ¥- 


X't fs-jsitfa 


(6) 


if* X- y*ts7S-jfa e&sJtfi) 

ZCx*~ z *fy 4 - %,Z + ?*j] h 
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The velocity induced by the serai -inf inite trailing vortex A«= is analogous 
to that of A'. Therefore; 




Hex*- 




? rf l7 (IIX- *- Jl 


+']£■?*&+&-*$ (9) 


If the Goethert transformation is made and the coordinates 
nondimens ionalized in terras of the serai -width of the skewed vortex, the 
following expressions for the nondimensional induced velocity ratios 
are obtained . 


4-7T ^^2 j \£, <frr J ^ (10) 


where 







■/- 





£=" " 




•r<n^2 

L 7? 2 


CTo+O CX 





( 11 ) 





g-'X^ 


+ 0 


+0 


*■%> 

K 
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and 


7?'* -fey* *?'l7 

^ p ' ■ *,>J 

+ ?'y 

/?y v 2W<^.;V cy+ jy + rU y * 

_ }f >*• y /■ ( / +■ <£*?-*s* 44- £ 'J 

2 v rs&ofe)* p ey+*S* s- p'I/^ 
y' / 

/vv ?-y* 

y^ yyyy^y 
■>"- C*'^ ?'!/ 4 


(12) 



y"^ y'rsss^fa — j 

/&f - )V<f y- //v 5 '7^ 


•=£ = 


X' — s^/ 


[C* f - &9nJ4t, t 'f‘ s (y '-'f- / 2 ' !/^ 


l 6 o 


. r i ' v 

The non-dimensional, coordinates (x , Y , Z ) are defined as follows 




■g ‘j-y- J 


7 = 


- >t 


*? - 


5 'V- 


( 13 ) 








The equations for the velocity ratios can he simplified for the following 
conditions . 

2 2 

1) if R2 or K 3 equal zero, then the terms in which these 
quantities appear can be neglected. 

2) if 0 £ .<:‘ = $<><: then the third and fourth terms of E^. and Ey 

can be replaced by 


7? 2 


/■rj 




and 


C ... ^ 

^ //V 

where 0 - 0^. = 0^ and R 2 = R 2 2 = R^ 2 


It should be noted that the condition 0 4 - 4 - - 0 Dt - is satisfied at 

every span station except where there is a discontinuity in the rate of 
change of leading or trailing sweep. 
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The equations can he specialized to that used in near planar lifting 
surface theory* In that case Z f is assumed small enough that only the upvash 
is considered in the determination of the wing vorticity. The influence 
function then becomes: 




'-/ llt*'- 


( y~o L J' /r 


J " v’** ilOc'+e****# 



/ yVy f- (/ 

(/' *5 ) lit* *+ + ty v /j 2 J ^ 


yV 7 

^ y'VV 


(i*0 


c?- V^mT) l [s z * r V* 


y' * jt't&jtfoi - c/*-t**~*&d 7 

+ Cy'-tF J H J 


The third and fourth terms in can be neglected when (x* - y 1 tan 0^- ) 
or (x' - y 1 tan 0 oi ) equal zero respectively- Also, when 0^> = 0 O . the * 

third and fourth terms can be replaced by: 1 


O'- *J LZjC/'t&a**?? ^yVx/7^ / 7j** m J 


where 0 
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The components of velocity induced by the skewed vortices can be written as 
the product of vortex influence matrices and the vortex strengths. 

{ 4 } - [ « li 1 

{5] « li] 

and 

^ ](-} 

where the elements of ^ j &y 3 and are computed from 


/Sr - 


( 18 ) 

tTTtfy. & 

*r * 


(19) 


and 



* 

11 


(20) 



(15) 

(16) 

(17) 
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Appendix D 


QUADRILATERAL VORTEX INFLUENCE EQUATIONS 


The Blot-Savart lav can be used to calculate the influence of a finite 
vortex segment on a point in three-dimensional space. The incremental 
change in induced velocity at a point in space due to an incremental change 
in length of a finite vortex is given by the following expression. 


dq. 


K COS 0 d 0 
4wh 


( 1 ) 


where; 

K = Vortex strength 

h = Perpendicular distance from the vortex segment 
to the point in space. 

0 = Angle betifeen the line formed by h and a line 

from field point to a point on the vortex segment. 

q, - Velocity induced by the finite vortex segment per- 
pendicular to the plane formed by h and the vortex 
segnent. 

A vector expression for c[ can be determined from figure (D-l) : 





POfNT OM VOFTEX 


Point 

tNFt.UENC.ED 


Figure D-l. Velocity Induced by Finite Vortex Segment 

1 6k 



She magnitude of^the velocity £ induced at ( j Vg ■> ^3 ) "hy the 
vortex segment s is given by the following equation after equation (l) 
has been integrated from jsf. to <?£■ 


If/ = o 


( 2 ) 


“where 


Ccrs-p 


<& ■ ^ 


CoS e < 


g 

/*///%/ 


She vector A is determined such that it satisfies the conditions of 
being perpendicular to s and equal to the vector sun tT « ^ - ct.§ 

■where "a" defines the length of h . 

Since; 


h 


- 7^ 


(3) 


and 


_!W, 

h ■ «■ 




then 


^'.g = ^S r — As r .s‘~ 0 (5) 

therefore; 


«. = 


• 5" 

r -«P 

g . s 


( 6 ) 
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After substituting “a" into equation (3)> ^ is defined as; 


h * 


— 7Z, f -S' 

w ~ 

X =■ • s 


(?) 


Also, a unit vector g in the direction of J? is seen to be equal to; 


XI __ * -s 

% x*/ 

The magnitude and direction of % are then expressed in terms of the 
coordinates of the control point ( x g j 5^ , ^ ) and the end points of 
the vortex segment { x ( t y t . J ? c . ) and ( ^ , X / ^ ) * if T, T 

rind K are defined as unit vectors in the x , y , and directions 
respectively, then; 



- (Xf-Xi)? J? 

*£ - ^ = tX'-Jtf)? 


and 


- Cx A -x f )i> 


( 10 ) 


The value of Tf a" is then expressed os; 


<z s 


fZf. • ? 

— * 
S • S' 


yaff ■ -X t - - K-H-ffo: 


(ii) 
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pnri the components of h by. 




*“ °W«:^ ?‘C +CV 2 tfjS* (12) 


fc, 


-/ t XX r ^J +*^V>->5XYJ: -y;-; +£^-3^^-^.; 

+*ty f -r t o’- 



The magnitude of f is found by substituting h *■/*/ - J<hf +<hy + 
and the following expressions for g°* e*. and c.*>* % into equation (2)* 


^.cc = i-p 


1 \(H - x ; f 7f *(r f -r t ft-? ! (*,-t c f i <*#-*? >S-?b f-Yf T*$%r %f 


( 13 ) 


The components of the vector f axe then given by the multiplication of 
the components of equation (8) by /£/ . 

o W '#-*X - C 3> X# - *;/<?’ 

Zit “ 7^ — — 

//2X£?/ 


(1U) 

/$**! 
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■where 


Zx*$Z ~[[CYf 


-ZtXXf-yjjJ \ 


+-Z Ct t-* 1 Xr,->i)-ir' L ] h 


(15) 


The velocity induced at a control point “by a vortex segment is then given 
by equation Since a curved vox*tex can be represented by a number of 

straight segments, this equation can be used to compute the induced flow 
produced by a vortex of arbitrary shape. 

The components ox velocity induced by a quadrilateral vortex can be 
written as ratios computed by the product of influence matrices and the 
vortex strengths. 


c* j • r 

* ]{i] 


( 16 ) 

feM 

+ Ml 


(17) 

and 





+ ]{t) 


( 18 ) 

where the elemei of X* 

i’ollowing equations. 

f /3y t and X-g , are computed 

from the' 

/* = 27 

■?rcei>7£- - % X& -Z) - d?/ * )J 

(19) 


? TT/hJjfcX^i 




(20) 


4-iT /h//7ZY^ / 

“ \ i 


168 



and 





? r ae ^iB - c - XfX >/ “ r<) -(Yf~ r £ )C^~^)J 

4 7T/Z// 75 


( 21 ) 


The£ sign indicates that the contrihutions from all of the sides of the 
quadrilateral vortex are summed. 
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Appendix t 


SOURCE FRUSTUM INFLUENCE EQUATIONS 


flow through nacelle Is represented "by a distribution of source 
frustums on the surface of the nacelle. The expressions for the velocity 
induced by a source frustum are given in reference 22. If the nacelle is 
defined by a series of points ( Yk ~ ) j 222(1 these points represent 

the ends of the frustums, then the midpoints of the frustums are given by; 


**i - 


XI £+ . -hXIt 


(1) 


YZl = Yi c „ + n L 


( 2 ) 


The surface lengths of the frustums are defined by; 


•= * t"i,i - ri^?] : 


(3) 


and the Xo component of the vector normal to the frustum surface is 
given by; 


" 




•=. — 


w 


2he component of the frustum normal vector is given by; 




xt £ +, - xt I 


COS’ oi^ L 


(5) 


If the points being influenced are denoted by CXz t J Y ? t -) , then the 
velocity ratios induced by a frustum of unit strength at (XZj t yz.) 
in the and Y*j directions are given by; 



“ - 4 





c/sr. 

j 


( 6 ) 
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and 



r 




CT) 


respectively. 

"Hhere, frcea reference 22.; 



/ 


[Cr2 ; + Y}) 1 +■ -X; 


( 8 ) 


and 


^-]cre^ * /+ (* V *j # * [ K( *‘ y) 

Also -where; 


* 


r*/- x- -o&t-sj? 1 

C*t t - r^i-Grs-Sj? J 


(9) 


**■=* Atj f'S; <2<*S« 


( 10 ) 


Yj - Vlj 4- s j S/aJX, 


(11) 


and Sj is the surface distance measured from C*fj ^Yij) along the 3 xh ’ 
frustum,. K Os*; ) and E C*,:,* ) are complete elliptic integrals of the 
first and second hind, respectively. Where, 



For the special case -where Y* t -) a&d rzy) are the same 

point; 


C%) e£ - +- ., ^(.(3 <-4|- 


(13) 
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and 


(.—) - 
' vi Ju 


- zrr cos =<. 


v ->-Jn §' ) s .- 


■C 3 * ^ -g 


where s t - - £ -OS 


for « £ >.o 8 , set --^6 and use; 


(f — 7— J., — “2ll CtJSflifj' h ^ +■ JZn -i; ) 

% \c» 'a 1 0 


~(3 + 3A-,§ - 3^. 2 s^.)S 3 


r-as,: 

- 2 Uv^’ -zl < 


f . . «*v 

£ i*4 4 


(-£ ) t .. ~ + ^^.cos^. [zs. C-? * ^5* 

"U - 'I ^ 


■where 


ds = *Vz 2 -^^6 


<1^ 

—°8rz c 


- + 'GZYZi 



The velocity ratios in the Xjj and Y^j directions are then given by: 


M - , 

r« mi 

(IT) 

and 

&H 

> Is] 

( 18 ) 

■where the elements of -5^ and are] 



( 19 ) 



( 20 ) 



APPENDIX F 


WING AND °yLON VORTICITX SERIES 


A modified double Fourier series suitable for representing the loading 
component of vorticity over a thin wing in subsonic flow will be derived in 
this section. The series will be general enough to represent the loading 
over a wing of arbitrary planform, twist and camber, and with leading and 
trailing edge flaps. The boundary conditions of zero net load at the wing 
tips and trailing edge are satisfied by all the terns in the series, 

A series applicable to wings with continuous spanwise variations of 
leaair.g and trailing edge sweeps and chordwise variations of mean camber line 
slopes was first developed by Blenk (l2 ) • The series was derived by combining 
tbe expression for the spanwise variation of circulation developed by Betz (13 ) 
with the expression for the chordwise variation of vorticity developed by 
Birnbaur. (l4). The expressions by Betz and BIrnbaur. were derived by means of 
thin wing section and lifting line theory, respectively. 

The general expression derived here will be developed In the same way 
as that by Bleak, An expression such as Bleak's will be derived first and 
then modified to represent planx’orms of arbitrary shape and wings with 
leading and trailing edge flaps. Even though the expression will be derived 
for a thin wing in a uniform flow it will be applicable to a wing-pylon -fanpod - 
nacelle combination If the induced flow from the fanpod and nacelle is smooth 
and continuous and if the coefficients in the series are determined such that 
the boundary conditions over the wing, pylon, fanpod, and nacelle are 
satisfied simultaneously. 

In thin airfoil theory, where the camber is assumed small enough to 
replace the singularities on the mean camber surface by singularities in a 
reference plane, the velocity Induced by a chordwise distribution of vorticity 
k' (x) is given by the following expression. 


/ 


v t - C x t, 1 - 2 mr 



X'C/) dx 



( 1 ) 


Where x is measured from the leading edge and x Q is the location of the point 
being Influenced. To satisfy the boundary condition of zero net flow tlirough 
the mean camber surface of the airfoil, the sum of the normal components of 
velocity due to the freestream and the distribution of vorticity is set equal 
to zero. The following constraint equation results. 



( 3 ) 


where a is the airfoil angle of attack and st is the slope of the mean 
camber surface. Substituting (l) into (2) gives; 


sir 



YU) <** 

X D - * 




o'* J 


a 


This equation can be transformed into polar coordinates using; 

X - \0~ Co-Sdp) 
to obtain 



Cos & - '<s> 0 


VS« C ex 


c/g 




and solved for . ifte flat plate solution is already known from the 
Joukowski transformation of the flow over a rotating circular cylinder such 
that the Kutta condition is satisfied. This solution is given by; 


= -2 eor<sy z ( 5 ) 

where is a constant . It can be shown that #0 not only equals ex for a 
flat plate but also equals <x in the case of an airfoil with camber which 
is symmetrical about x * | • 

in general the distribution of vorticity over a cambered airfoil is a 
combination of a term such as in equation ( 5 ) and an infinite series which 
can account for any arbitrary deviation in vorticity about the ccr^z. term. 
A series which satisfies the Kutta condition and when combined with the 
Cor<s>/-z. term is sufficient to represent the vorticity over any cambered 
airfoil is the following infinite series. 


Y,CS0 » -2 v£, 2^ ^ 


(6) 




If equation (6) is substituted into equation (l) 


r 77 * -^><p . v/ . 

s*i-~ dO Is. (_ & \ 

zr J dc-3 (9 _ COS <3?, V -~ 

4? 


(T) 
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( 8 ) 


and integrated the following expression for 



is obtained. 


% 

v„ 


f . y — — 

" Z 


sJ~, 


43^ C os’ <9 


Since equation ( 6 ) and equation ( 8 ) are Fourier conjugates, equation (7) can 
be inverted to obtain a general integral expression for Y z (s>) in terms of 
any arbitrary ¥* cs> a ) which satisfies equation ( 8 ). 





4- 

5^3 

CCS <S> — CC>S <S> a 


d& 0 


(9) 


This expression cen be used to reduce the requirement of an infinite 
series, to represent the vorticity over airfoils with discontinuous free- 
stream slopes such as in the case of leading or trailing edge flaps, to 
that of a finite series. This is done by representing the net pressures 
due to a smooth camber and angle of attack by a combination of and a 

finite form of Y z os>) such that; 


Y&) - S^fcso + ^ 09 ) - /- 




V-t 


(10) 


then adding to this expression additional terms evaluated from equation ( 9 ) 
for distributions associated with leading and trailing edge flaps. It should 
be noted that both and ) are needed to represent smooth cambers 

which are nonsymmetric about x- , even if equals zero. This is because 
the antisymmetric sine terms in equation ( 10 ) contribute to the value of /?*,. 
For the same reason both y r d®j and along with the additional terms 

from equation ( 9 ) must in general be used to represent either leading or 
trailing edge flap deflections. 

The relation between the sine terms, ci , and /* 0 is derived in reference 
end repeated below 






/ 


4sm 7 - / 


(XI) 


This relation demonstrates the dependence between the terms in equation (10). 
This dependence could cause the Influence matrix, eventually to be inverted 
to solve for the coefficients /9 e and } to be^illconditioned . In the 
limit as « it can be shown that cor<o/ z =■ ^ e? for c><&< rr 

which would make the influence matrix singular. This can be avoided if the 
influence is considered over the complete range <os & z rr~ as done by 
Wagner ( 15 ). 


Due to the finite nature of the series in equation (10) and the con- 
venience to he derived from its use, equation (10) and its modification to 
account for flaps will he used here . It is thought that the above mentioned 
problems can be minimised by judicious selection, of the points where the 
boundary conditions are to be satisfied. 

The special chord’, rise terms needed to account for trailing edge flaps 
can be derived by setting QC<9 0 J - o for and ^ for 

tys cSi s xt~ in equation (9) and integrating over the chord. \lh ere <s>+ is the 
trailing edge flap hinge location and g the deflection 


=r 



S/ 


<d <S > 0 


( 12 ) 


7T 


r <- C<S> S <£ 





rr 


( 13 ) 


A similar expression can be obtained for leading edge flaps by setting 
•*. - ^ for o< <s> 0 s and = o for <S‘ K <s>„s rr } where <s>^ is 

&e leading edge flap hinge location and £*. the deflection. After 
evaluating equation ( 9 ) the folloirilng expression is obtained. 




YT' (£-/sO <sw / 


( 14 .) 


The subscript k is used to indicate a leading edge flap or droop, both 
can equally well be represented by equation (l4). 

All the expressions derived thus far have been for a tvo-dimensional 
airfoil. A set of representative spanwise functions, to represent the span' 
vise variation of circulation over a lifting line, will be derived at this 
point. These functions will then be combined with the chordvise tenr.r. to 
represent the variation of the loading component of vorticity over a finite 
lifting surface. 
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In lifting line theory the induced velocity due to the spamrf.sc vari- 
ation of circulation is given "by the relieving egression. 

~T^vr~ *' C 1 ?) 

- d / 2 

This expression can "be transfer. :ed into polar coordinates by letting 
j where ^ ■ 


. L C ir C* P /c*4 t ) 

*=■ 2 TT^> Cos jzf - CoSf* 


d £ 


(16) 




A series which satisfied the boundary condition of zero circulation 
the ving tips and is sufficiently general to represent any spamri.se 

-.j _ • .j. • ■* *,^* **.- 


variation of circulation is the folio - , ring sine series 


^ 


oQ 

Z. 

uu- / 


‘S^ 0 S’z/v' 


(IT) 


If equation (l6) is multiplied through by sin $ and equation (17 ) 
substituted in for the folio - , ring convenient integral equation is 
obtained. 


—r <.$) 



7T tip, 

C.o-5 _ 


■r - /*£$) y 

c~ 2*0, b J 
G&sS ^ 




(18) 


since : 



<30 

y yJ S^J C.o^ wv £ 
Wc/ 


(19) 


and that after substituting equation (l$>) into equation (18) and integrating 
it is observed that tj and Cs*J.5*w<^ are Fourier conjugates. 


*o 

*= 7F J CDSxfy - <LC>*S. 


s y~ oo 3^ *5 “ ‘/■'J w 


( 20 ) 


"»s/ 
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(21) 


It is also observed from equation (20) that 


OQ 

£<•>.)= IT 

r -* *)~i 


V ^> gju 'SssV V*J 4> 
~5"/sXJ 


Co 4- £ ^ *** 


which demonstrates that equation (17) is sufficient to represent any span vise 
variation of twist. 

Since J^C^I^vj^and gxe p our i e r conjugates, equation (18) 

can be inverted and a general integral equation obtained for calculating the 
spanvise variation of circulation due to discontinuous rates of change of 
induced velocity such as obtained when the spanvise variation of leading or 
trailing edge sweep is discontinuous or partial span flaps are deflected on 
the leading or trailing edge. When these special spanvise expressions for 
circulation are added to equation (17), the infinite series in that equation 
can be reduced to a finite number of terms . 

When equation (l5) is inverted the folloi/ing integral equation is 
obtained . 




(22) 


This equation is then integrated with respect to $ to obtain a general 
relation for calculating the spanvise variation of circulation due to any 
arbitrary induced dovnvash distribution. 


ft ,/) ' (' 7r ^. ftp ) ^ 


(23) 


The type of dovnvash distribution due to leading and trailing edge flaps 
is Ino’.m, however, it is not obvious what type of downwash is produced by 
breaks in the leading or trailing edge . This can be determined by computing 
the variation in dotniwash, over the planfom of interest, produced by a vor- 
ticity distribution vliich •\rould produce a constant dovnvash over a high 
aspect ratio elliptic wing. Any variation in dcrc/nvash, relative to the 
constant distribution on on elliptic wing, can be associated with the dif- 
ference in planfona shape between the wing being analyzed and the elliptic 
wing. In this way the extent and shape of the dovnvash discontinuities due 
to the breaks in the leading and trailing edge can be studied. It is only 
necessary to account for the discontinuities in the dovnvash distributions 
by means of the special spanvise functions because the smooth variations 
in the downwash will be accounted for by equation (IT ) * 
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The special spam/ise functions necessary to account for these discon- 
tinuities are then variations of circulation which will produce downvash 
distributions with the same shape as the discontinuous variations found 
using the above procedure. It has been determined that if the discontinuous 
deviations are represented by polygonal dorm.wr.sh distributions, that the vari- 
ations of circulation due to these downvash distribution calculated freer 
equation ( 23 ) , do an acceptable job of accounting for breaks in the leading 
or trailing edge. Even though the discontinuities are of a logarithmic 
type, the polygonal downvash distributions are acceptable due to the presence 
of the other terras in equation (l?)* Also, it is felt that the use of the 
more simple polygonal functions is justifiable because of the eventual 
representation of the vorticity by a finite number of discrete vortices. 

[These polygonal do van rash functions can be built up by superimposing two 
basic types of distributions. These basic distributions produce variation of 
circulation which rail be called "M n or IkuLthopp functions, named after 
Hulthopp, who first derived the functions, and "p" functions because of the 
polygonal nature of the dormvash produced by them. The special functions 
necessary to account for flaps, leading edge droops, and ailerons can also 
be built up by different combinations of the Kulthopp functions. 

The Kultliopp functions are associated rath air induced downvash ratio 
) which is zero for -/£ ^ and unity for - S • Substitutirg 

this dovnvash distribution into equation ( 23 ) and integrating by parts the 
following integral relation is obtained. 






The integral; 


1 


£ 05 ^ — ' 


i dL^ric4, - '"Tie*, 


(ait) 


( 25 ) 


vhle'h when evaluated gives 


= cos a j + 

) 0 acs $~ Cc > s f ' 1 C{6*- / 


( 26 ) 


Therefore; 



^ & -s~/^ 


ISO 


( 27 ) 



Another type of I-iulthopp function necessary to "build up symmetric and 
antisymmetric loading functions is one in which the ratio ^ C^) is equal 
to one for and equal to zero for *%*< -? < / . For this downvash 

distribution the spanvise variation of circulation is best calculated by 
transforming the coordinates such that; 


$ - - /0O 


( 28 ) 


Then equation (23) becomes; 

.ref) \ J_( ** &'"*$.<■£'- 12- l ^ A * 


(29) 


~ rrj j X evisf'-a&'Sft' 


(30) 


Since equation (30) is identical in form to equation (24 )j the following 
result can be obtained in the same way ns equation (27)* 







~> 


(31) 


Then after substituting equation (28) into equation (31) the following 
expression for the port wing panel is obtained. 


\2bA. 


~r) - j- fecc^S * gP5^i ) ^ l e jf- y 


(32) 


A symmetrical and an antisymmetries! I-Iulthopp funecion can be formed 
by superimposing equation (27) and equation (32). The symmetrical function 
will be used when the wing contour is symmetrical about the plane of pitch 
and the airplane is rotated in this plane. The antisymmetric function is 
lined ifhen the ailerons are deflected and the airplane is in a rolling 
maneuver. The symmetric Multhopp function is; 
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rftfj <*% = £[ cct=>s ***& 


C$*- <zO 






(33) 


The antisymmetric Multhopp function is as follows. 


, f is>^^Cf>’ , ~i s i / 

^*0* ’ £/<“*/*- '/ 


_ <r<?^r ^ v* <?*=>; 




( 34 ) 


The "P" function is associated with a downvash distribution which is zei’o 
for '/<■ ^ s' ^ ^ and varies linearly from zero to unity for -£■*<•<£ s*/ - 
For this distribution equation ( 23 ) becomes; 



/ — fie-s - 






e/6 s 

i* £>J / 


( 35 ) 


After integrating by parts; 


' rrO- &&&$*) C 




/rs?^v ^ 


'^s^tjCfS 




- -;^s> 


( 36 ) 
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The integral; 



— CL ~Z, 




^ 8 


r &* 

j f e-c^ ceT iCf. -'£) 



( 37 ) 


whexe. 


j? f-y '$) f 

j >S z<$ *&rf C4-& °# - ^ e-°T-£C$+<t^ - i cos 

_ - Z $*s>*3zf> (38) 


therefoi'e; 


(**L) = 


/ 


r~ 


277 cs- ao&fS^j 


l c 


C^c>S <pS** <ZC, 


*<*■> 1 */ 3 ZJW&* 7 df 


(39) 


+ Cz ff*- S>s*f& - ■£ ff*sr/sO 2 j^? 


Then similarly the "p ,T function associated vith the port iring panel is 
derived "by making the sane coordinate transformation as used in the 
Multhopp function. 

Let = 180-^, then; 


,£<£L\ _ ../ (kO) 

<^*,4r “7T J C 
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Since equation (40) is identical in fora to equation ( 35 ) its integral is of 
the same fora as equation ( 39 ). 


27ro-<sos<f) 




7<p*C<>S ^ *- S ' jW 0'-* *4* J (4l) 


After substituting 0 - 180 - 0 into equation (4l) the follo’.'ing expression 

is obtained 


2 £»*£#■ /jp 


^ ,7- ,„s- fC C &*+■<?) f 


•znCt-ecsj+l 
r 1 ( 7 . 


(te) 


As in the case of the Multhopp functions a symmetric and an antisymmetric 
"P" function can be formed from equation ( 39 ) and equation (42). 

Fne symmetric "p" function is formed by adding equations ( 39 ) and (42). 

PC6<*%= “ JfFc^ct^^y 


r jfc 


f-(_4 ^ S' 


,7 


m 



The antisymmetric "P" function is then obtained by subtracting equation (42) 
from equation (39)* 


#*) ■=. li cos d*- ggg l e&e ( 

-zrrc/- ce&ft*) (J “ Aw £ <■ &*+ ftJj 


, £ 2 [COSftf* 

- s-"^ z tfij- 

The spcnwise variation of vortlcity necessary to represent polygonal 
downvash distributions are then built up by the superposition of th gs?C(% d") 
and "Ft#, ft 1 *) functions. Tlie proportion of each used depends on the location 
of the break in the lifting line ^ and the ranges lb and over which the 
downvash is affected. Tile range varies vith the degree of discontinuity in 
the sweep angle of each lifting line. In general; the change in sweep angle 
is not the same on the leading and trailing edges, in which case P^ and Rq 
becomes functions of x /a . However, the special functions are calculated for 
constant values of Pi and P Q with the assumption that any variation iriLth x/c 
can bo compensated for by the other terns in equation (17)- This appears 
to be a safe assumption since the range has been observed to vary only slightly 
about the value .05^- The actual plot of the induced downvash can be obtained 
from the elements of the influence matrix. 

The special functions which must be added to equation (17) to account 
for discontinuities in the variation of leading and trailing edge sweep are 
given by the folloiriLng relations. 



For ^2 - G 


■pc#) =, WJ - 


4 -pc^W) * 


ispL ^s"/^ ; 


(^ 5 ) 



For ^ £ - 


r- £ rtf*, 

- PC j ^ ) 
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f " A I / “ K " ^ ces't^*# 
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For J %£/<%,/* 0-*„) 


-pet) - ac&'c^-^)) 


~ PC *' 40 


-^o )C/7^)7>C^ G&~? ^ /- )) 


W 


For Cs-Z t ,)$/'%'/<r/ 


-pc<?) ~ Q- %, s&iJCM') 7?c ^ 


* <' / -^;/^;. J/ 7=^ <&. ; 


(*8) 


These sane relations are used for “both symmetric and antisymmetric loading . 
The only difference "being the use of f^C^j <^*) s and 7 1<&; for the 

symmetric loading and 7^ for antisymmetric loading. 

Also, equation (45) is always zero for the antisymmetric loading. 

The special spanvri.se functions for flaps, leading edge droop and 
ailerons are similarly developed. The function added to equation (IT), in 
addition to the above relations, for flaps and leading edge droop is as 
follows : 


« /*7C X> y ^ ^ ) 


m 
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fa and are the inboard and outboard spanwise locations where the partial 
span flap begins and ends, respectively- For ailerons, the special span- 
wise function is the sene as equation (49) except that the antisymmetric 
Eulthopp functions are used instead of the symmetric. One special function 
Should be used for each flap or leading edge droop for which and ^ are 
different. 

At this point all the necessary expressions needed to represent both 
chordvise and spanwise variations in the loading component of vorticity have 
been derived. The expressions have been derived from lifting line and tiro 
dimensional thin airfoil theory. However, they can be combined to represent 
the vorticity over a lifting surface since their shape is representative of 
the type of vorticity expected over a lifting surface- The amount of each 
function is determined by satisfying the boundary condition of no flov through 
the mean camber surfaces of the ^ring at a serieG of points, called control 
points. In the case of a lifting surface , however, the velocity induced 
by the vorticity is computed using three dimensional influence functions 
instead of the two dimensional functions used in the derivation of the loading 
expressions. These three dimensional influence functions are derived in 
Appendix C. 

The combined general expression for the vorticity over a thin lifting 
surface is then obtained by merging equations (10), (13), and (l4) with 
equation (17) and the appropriate versions of equations (4^), (46), (47), 

(45), or (49). After transforming the spanwise coordinate bach to ^ ^ lG 
following expression is obtained. 
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■where, for a symmetric span load; 
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and where /^ equals the total number of span loads, ^ equals the number of 
special span functions needed to account for breaks in the leading or trailing 
edge and for flaps and leading edge droop and where the ’%> { - < Z)' > S are the 
special spanwise functions. The coefficients , and ^*^>are 

the unknouai ' s to be computed such that the boundary conditionr arc satisfied. 

If the loading on the irlng is antisymmetric change the poorer on ^ 
from ) to -z^“/ and use antisymetric?far,cti one. 

Since the pylon does not utilize a spanwise series or have flaps, the 
vorticity is defined by a chordwise series only of the following type at 
each spaa station. 
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rn tnis case F 0 and Fjj are the unknown coefficients to be determine by the 
■-.olution of the influence equations. 
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WING AKD PY ! OI’l VORTEX STKi.'J'TTHS 


The expression for the vorticity over a lifting surface Is given in 
Appendix P oy equation (J?0). The unknown coefficients of the vorticity series 
are determined such that the boundary conditions of no flow through the 
surface of the complete configuration are satisfied. This is done by setting 
up and solving a set of simultaneous aerodynamic influence equations which 
represent these boundary conditions. In order to set up these influence 
equations it is necessary to sum up the velocity induced at a series of 
points by the vorticity. This summing process is in general the integration 
of the product of the vorticity and an influence function over the entire 
lifting surface for each of the boundary points. 


In order to eliminate this analytical integration the continuous 
vorticity distribution is replaced by a grid of constant strength vortices 
as shown in figure 4. This is in essence a numerical integration scheme 
which does an excellent job of representing the induced velocities with a 
minimum of numerical computation. The technique of representing the vorticity 
over a lifting surface by constant strength vortices was used by Falkner in 
reference l6 . Falkner utilised horseshoe vortices and obtained good results 
for wings of low or moderate sweep. However, the number of vortices required 
increases for wings of higher sweep and eventually the representation 
becomes unacceptable. This is because the bound segment of the horseshoe 
vortex is perpendicular to the free stream, rather than along a line of constant 
vorticity. Since lines of constant vorticity tend to lie along lines of con- 
stant percent chord, the so-called bound portion of the horseshoe vortex has 
been skewed to lie along constant percent chord lines in this analysis. This 
procedure should also produce a better representation of broken leading and 
trailing edge planforms. 

It can be shown that discrete vortices produce the same downy ash, at 
points half way between the vortices, as a continuous distribution of vorticity 
represented by small linearly varying distributions over regions equal In 
length to the distance between the discrete vortices. Since any distribution 
of vorticity can be represented by a set of linearly varying regions, the 
number of which depends on the rate of change of the gradient along the 
distribution, it is reasonable to use discrete vortices to compute the 
induced velocity produced by a continuous variation of vorticity. 

If the arbitrary vorticity distribution, as shown in figure G-l , is 
represented by both linearly varying vorticity distributions and by discrete 
vortices the following expressions for dovnwash at a boundary point midway 
between the vortices are obtained. 
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Figure 0-1. Vorticity Distribution 

She downvash induced at x = 0 by tvo dimensional discrete vortices at 
x = - 4. and x = equal in strength to /J* and , respectively, is 

given by: 
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If the continuous vorticity over the intervals -/=S * S ° and 
is represented by Taylor series expansions about x = - and x = 
respectively, then 

* c?) =. w(-4) + <•* + b y'c-4J + 

aw = vck) *■ cx- 4; ^(4) +■ • • * ( 2 ) 


If these series are tnmcated after the first two terms, the vorticity is 
represented by tvo linear functions whicli produce the following do^mwash 
at x = 0, 
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Therefore, after solving equation (a) for ¥*£{) and y'ci) at :c = 0 and 
substituting into equation (if); 
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Then in taking the limit, it is observed that; 




Therefore; from equations ( 5 ) and (6); 




which is the same result obtained from the discrete vortices. 


( 8 ) 

C9) 


Even though this argument and others to follow in this section are 
based on the two dimensional specialized Biot -3a vert relations it can he 
shorn that there is no loss of generality and that the above conclusion and 
■the relations "between the chord vise variation of discrete vortex strength 
and the continuous distribution of vorticity, represented by those vortices, 
are valid for finite vortex segments and vortex grids such as shown in 
figure 4. More general expressions which give the induced velocity due 
to skewed vortices of finite length can he derived from dowuvush relations 
given in Appendix C- Using these relations the lo*., awash produced by the 
bound portions of the discrete skewed vortices, vx a point mid -way between 
them, is given by; 
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Where S is half the y distance between the ends of the bound vortices. 
The streep angles of the port and starboard portions of the bound vortices 
are represented by 4 and. do , respectively. 


For a continuous distribution of vorticity the downvash at x - 0 is 
given by the principle value of the following integral. 
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where the variation in sweep of the vorticity is represented, by a Taylor 
series expansion about x = 0 . which is truncated after the first two terms . 
Therefore j 
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where 
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After substituting equations (lk) and ( 15 ) into equations (11) and (12) 
the following influence functions are obtained. 
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where 7? -=z t=s%* + j ^ ^ ^ ce>) ^ ;T =&£&&£&) 
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If equations (l6) and ( 17 ) are substituted into equation (13) the 
following integrals are obtained. 




f(Zx+*7 a ? + 


OTo 




* £ 




7c/x _ £ 

i+sKF+eV 4/r ^° h ^ x 


- r/ J^;. 




lWr- x : 




. 7> 


r°r^ -, , ')C%z'x+Cz%+*>%‘ xz'x +cz%~/)t„ 




4- 


& 


*Jc? c x **■ y -Jc%x ^ t p) 


/^£> ■y'. "5 

• ~' T^ ~; - — - i- 7k. - 7 : J-T'H* 

^X-hsT^s 1 - X X c »] 




^ttKo J 0 (- / ? J C^xV^V ss* J<5%x + ST 0 ) r i- s* 


y- 




-7- y 






+ IL - ~r~o 


-h v‘ + n'J cJx 


( 18 ) 


where means to take the principle value of the integral. 
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The integrals in equation (l8) are evaluated as follows; 
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and 

-£ j, 

^ * c * )T °tzrS it- rco ^ £z>[ % = ***>£ fxi *°**/r/ 

go logarithmically to infinity as e-^o , This is a consequence of eval- 
uating the downwash at a break in a line of constant vorticity. If the more 
coi:.plete planar influence equation, equation (ll) of Appendix C, is used 
in equation (10) an additional term is obtained in the limit as e -♦ o. 

This limit is difficult to obtain because o? ^’ 4 - and , however, it was 
noted in reference l.j' that if the downwash were computed at a point on the 
surface of the airfoil instead of in the X = 0 plane, that the downwash at 
Y - 0 due to the break in the constant vorticity line would be approximately 
proportional to the vorticity at that point. It was further noted that the 
downwash was not sensitive to the thickness ratio of the airfoil section. 

A more complete study of the downvash at breaks in the vorticity lines will 
be made, using equation (ll) of Appendix C, at some tine in the future. The 
results noted in reference 17 are not obvious when equation (ll) is used 
due to the complex nature of the integrals involved. These integrals may 
have to be evaluated numerically as in reference 17* 

From a physical point of view there is no reason to assume that the lines 
of constant vorticity will be discontinuous at any point on the wing except 
possibly at the leading edge of span stations where the sweep is discontinuous. 
It is therefore a result of not knowing the shape of the constant vorticity 
lines a priori, and the necessity of assuming shapes before they can be deter- 
mined, which results in the singular integrals 115 and I 0 -j. The actual shape 
of the vortex lines used becomes less important as the size of the vortex 
grid decreases and there is no reason to assume that a set of pressures, for 
which lines of constant vorticity are aligned with the assumed vortex lines, 
cannot be converged to after a certain number of iterations. For these 
reasons the exact behavior of the integrals 115 and I05 in the limit, as 
(X,Y,Z) — (0,0,0) is only of academic interest at points other than the 
leading edge. 

The complete integral in equation (l8) will be evaluated here except 
for I<5 and Iqc, which will be replaced by the values determined numerically 
in reference 17 * However, the comparison between dovnwash computed at 
(0,0,0) by integrating the continuous vorticity times the continuous influence 
function and that obtained by numerically summing the downwash induced by 
equivalent discrete vortices will be limited to the case where the lines of 
constant vorticity are continuous. For this case jzl — A ^ & and 
therefore 15,5 + I05 - 0. 

'Tow continuing with the evaluation of the integrals of equation (2.V-. 
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’When the integrals associated with the downvash due to the vorticity on 
the inboard side of the continuous vorticity patch are substituted back 
into equation ( 18 ) the following downwnsh express it a is obtained. 
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Where I^cj has been replaced with 77* ^ ^ from reference 17 ■ 
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An expression similar to equation (19 ) can "be given for the dovnwash 
due to the starboard side. 
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( 20 ) 


Equations ( 19 ) and (20) when added represent the total dovnwash at 
x = 0 due to a distribution pf vorticity, where the sweep of the constant 
vorticity lines vary with x. The span of the vorticity lines is arbitrary 
in length and can be vaijied to determine the effect of grid aspect ratio on the 
ability of discrete vortices to represent the dovnwash of continuous vorticity 
distributions. Also, the sweep of the port and starboard portions of the 
vorticity patch chn be varied independently to represent a kink or crank in 
the wing planform. Therefore, the ability of discrete vortices to represent 
breaks in the leading and trailing edges can also be checked with these 
expressions. 
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In order to obtain the effects of grid aspect ratio and vortex sweep 
on the accuracy of the induced downi-rash calculation, when the continuous 
vorticity is represented by a grid of discrete shewed vortices, equations 
(ll), (12), (19),* end (20) will be specialized to the case where 
T 4 ' = T 0 - 0, T d - - T 0 = T, and the vortex grid aspect ratio is finite 

The dovnvash at x = 0 due to two finite skewed vortices at x = -Jr- and 
at x - •— is obtained fron equations (10), (ll), and (12). 
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Where J5* - sec. ^ y f - and Ki and Kg are the strengths of 

the discrete vortices at x =.-•#? and x ^ ^ respectively. 


Figure ©-2. Skewed Vorticity 
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If K]_ and Kg are the absolute values of K;j_ and Kg, respectively, -where 




Then; 
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Also, from equations (13), (l9),j and (20) the downwash at x = 0 due to the 
continuous distribution of vorticity tf(x) from x = --f to x = ^ is 
given by; 
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It is interesting to note that the dovnwash induced at x = 0 due to 
two infinite shewed vortices, one at x = and the other at x = #z , 

and that induced at x = 0 due to a continuous distribution of infinite 
skewed vorticity are equal. This is observed in the limit as f-^<D 
in equations (23) and (24). 

For f - 0 the dovnwash ratio at x = 0 due to equation (23) is. 
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and due to equation (24) 
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whieh is the some as equation (25). It is also noted that for 0 - 0 "both 
equations (25) and (28) reduce to that obtained for the two dimensional 
case in equations (l) and (9)* 

Another interesting observation which can be made from equations (19) 
and (20 ) is that the downvash is independent of the vorticity, at the down- 
wash point, when the lines of constant vorticity are continuous. For the case 
in which the lines of constant vorticity are not continuous, the additional 
downwash is seen to be proportional to the vorticity at the downwash point. 

The ratio of equation (23) to equation (24) is shown in figure G-3 as 
a function of grid aspect ratio for two sweeps . It is evident from this 
figure that the simpler infinite aspect ratio expressions can be used to 
determine the strengths of the discrete vortices used to represent the 
continuous distribution of vorticity. These vortex strengths -will be 
determined such that the same dorawash is obtained at points midway between 
each vortex and the same circulation developed at each span station, when the 
discrete vortices are used, as would be obtained by analytically integrating 
the influence of the continuous vorticity. 
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Grid Aspect Ratio 

Figure Gr— 3 • Comparison of Downvash due to Discrete Vortices and 

Continuous Vorticity as a Function of Grid Aspect Ratio 



If S 
influence 


Is allowed to approach infinity in equations ( 11 ) and ( 12 ). 
functions for an infinite aspect ratio grid are obtained. 
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Ihe downwash at Xj due to a vortex K* at is then given by; 
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Ihe downwash at Xj due to the distribution of vorticity ytx ) between 
X.-£ and is similarly given by; 
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Ihe vortex grid is show schematically in figure G-I 4 . 
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Piguie 3 -4* Schematic of Vortex Grid 
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Equations (31) and. (32) are then summed over the chord to obtain the 
total downwash at each ,j point* The discrete vortex pattern produces the 
following total downwash at any j point. 
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A s imi lar expression can he obtained for the total downwash at any 
j point due to the continuous vorticity distribution. 
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If the sweep of the lines of constant vorticity are held constant 
over each of the grid panels 7 then; 
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Equation (33) can "be put into a similar form as follows. 
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The vorticity V (x) is defined in equation (50) of Appendix A. 
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where <cj — <?£^s''^/- ^ 


Therefore; ^ (x) is a lineal' combination of /t£ load shapes. 
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There will be one set of K L - T s for each of these load shapes, therefore; 


A t - - Kz + *- *> 


K-t. 


Si,-?. 


(39) 




209 


It is these values of /*■ which must he determined such that the 

downwash and section circulation are the sane for both the discrete and 
continuous loadings. Equations ( 35 ) and (36) will supply Mt 
equations. An additional Aic equations can he developed from the condition 
that the section circulations due to the discrete and continuous loadings 
agree. 

Hie total section circulation per freestream Wlocity for the discrete 
vortices is; 

= Z Z w 

£' = / 


The coirespouding circulation ratio for the continuous loadings is; 
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where <$* and are the leading and trailing edge flap hinge locations 
respectively. 

After substituting equation (38) into equation (35) and equation (39) 
into equation (36) the following equations are obtained respectively. 
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where 
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Since the elements of the matrix T§ t - 1 are independent, equation (42) 
and equation (43) are equal at the j points provided £ Qc^ Jf,-)- 
for of the chord wise load shapes. 

Also, since the functions E CX ij % ) obtained for each value of i are 
independent, then gTOtj j^)- S'C*,- ) for ^ chord ' ri ‘ se load shapes if; 



The matrix equation (44) represents *£ X*tc. unknowns and 
independent equations. An additional Ak equations are obtained from 
equations (39)> (4o), and (4l). 


['.TV] 


r AlO; 


I J S *rt ; "b ‘■'bo.-yjc J j 


t<~. 


K, 




= ° 3 








1 


/* Mr 




212 


are the a evaluated 
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since; 
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since; 
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After substituting equations (46), (47), (^)> (49)> (50)> and (51) 
into equations (44) and (45) the following set of simultaneous equations is 
obtained from which the discrete vortex strengths }c ct ' , Ki/S , /*/£ > 
and J^urt can he detenoined. 


~K,A 4, 3 3 * ' *, Xvw/xa* J *S'/uL, } ^'/u^ 
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t j -4002^ j *■ *.5 -eosif/<i-3^^- j — ^yr-a>c) 5 Or-<s> /c . ) 

hs n > ° ,••'••, 0 

I-Iote, in equation (52) the W,- / column is -^^5 for os^-cip- and <s>/ 
for <s»t < <s>j-s. rr . 


Also, the *Vu column is -<fr-dz*) for &<z <; <sj, c and sero for *£}•£■ /7'- 

in equation (52) 
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Equation (52) is solved for the ^ 

sides by; 


by p re multiplying both 


The discrete vortex strengths ore then given by; 
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( 53 ) 


Equation ( 53 ) is for symmetric span loadings, if the span load is 
antisymmetric the exponent zcy^-i) is replaced by ?w-t and the special 
span functions are computed such that they are antisymmetric. 

Since the pylon does not utilize a spanwlse vorticity series and doeb 
not have flaps, the discrete pylon vortex strengths are given by; 


= <‘-^^[c k '-‘/ uu )4 *■ zz c^u .) « 7 ( 5 >») 

^ ■nS’G.f 

vhere F 0 and Fpj become the unknown coefficients to be determined from the 
solution of the influence equations. 
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Appendix H 


FUSELAGE AND ^ANPOD VORTEX STRENGTHS 

The perturbation velocity due to the fanpod is represented by a vortex 
grid on the fuselage surface. Tills grid is made up of a number of closed 
circuits of constant vortex strength which are in turn represented by a 
number of constant strength vortex segments . The velocity induced by each 
of these constant strength vortices, at a point on or off the fanpod surface, 
is given by the algebraic sun of the velocity induced by each of the segments 
making up the closed circuit. The velocity induced by each of the segments 
is given by expressions derived in Appendix .D* The closed circuit can be 
of any shape, and can overlap in any arbitrary manner without changing the 
total induced velocity at a point due to the vorticity distribution being 
represented. It can be shown that there is a linear relationship between 
the influence matrices of twp different grid mappings. 

There is however, numerical advantages to one mapping over another. 

It is for this reason that a quadrilateral vortex grid with the smallest 
possible grid size is chosen. It can be shown that the Influence of a 
quadrilateral vortex tends to dump out at a rate proportional to ( ^ "f" 
where a is the vortex dimension parallel to r and r- ° /sV is the 
distance between the point being influenced and the centroid of the vortex. 
Therefore; the smaller A is, the smaller the area being influenced by each 
vortex and the smaller the 'number of vortices which need be considered as 
contributing a significant amount of induced flow at each influence point. 

This fact can be used to reduce the computing time per case by limiting the 
vortices for which a contribution is computed to those for which ( } z is 
larger than seme prescribed value. The vortices for which ( )" is 

smaller than that value should contribute such a snail amount that no 
discernible change in the influence matrix is made, if their contribution is 
neglected. It is - therefore advantageous to use a grid pattern for which 
the grid elements are nearly square and as small as possible. This is 
obviously not the case for the standard horseshoe vortex which trails to 
infinity behind either the wing, pylon or fanpod . 

Another method for reducing the computer time per case is to constrain 
the strengths of the quadrilateral vortices such that the amount of each 
constraint function, rather than the strength of each vortex, becomes the 
unknowns in the problem. In this way the number of unknowns can be drastically 
reduced, thereby reducing computer time. The constraint functions, which in 
this case are vorticity distributions, must be representative of the type of 
spatial variation expected of the vorticity over the fanpod, and must be a 
series for which end conditions, such as those at the nose and the base, are 
satisfied. The series must also be such that as the number of terms are 
increased the boundary condition of no flow through the surface of the Wing 
and the fanpod is better satisfied, in a least square sense. 
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These vorticity functions can be continuous provided the surface 
gradients are continuous. It will be assumed at this point that the surface 
gradients along meridian lines of the fanpod to be analyzed will be either 
continuous or that special vorticity functions will be added as part of the 
program input to account for any surface discontinuities. Under these 
conditions the following modified Fourier series is sufficient to represent 
the longitudinal variation of vorticity over arbitrary, but smooth, fanpod 
surfaces. 












/• ] 
i v-Z(Nv F i-l)+i 


( 1 ) 


where 

do ^'( /- ? * e /c p ) 


This series was used because of its general nature and its ability to 
represent the vorticity over both blunt and sharp ended bodies. One set of 
,T F" values is determined for each of the <£> roll angles. That is, the 
vorticity is only constrained in the ^Longitudinal direction since the 
nature of the vorticity in the lateral direction is more uncertain due to 
the presence of the wing and pylon. 

The first and second terms represent the vorticity over a pointed nose 
and tail respectively. These terms are identical to that used to represent 
the vorticity over a flat plate at angle of attack. The sine and cosine 
series is sufficient to represent the variation of vorticity due to continuous 
longitudinal variations of cross sectional shape. The last term in equation (l) 
is a set of special functions which can be input to the program to account 
for variation in vorticity due to discontinuous variations of cross sectional 
shape or any other condition for which the standard terms are incapable of 
representing. The sine and cosine series should be sufficient; for most 
shapes provided enough terms are used. The special functions can, however, 
reduce the requirement for a large number of Fourier terms and thereby 
reduce computer time. These special terms can be computed from either 
slender body theory or Sfrcrn axi symmetric body theory such as used in the 
Dougl as -He unarm, program of reference 22 . 

Hot all the terms need be used for every body shape, however, if more 
terms are used than necessary the answer will be unaffected, since a least 
squared solution is obtained for the influence equations. Computing time, 
however, can be saved by judiciously picJdlng only those terms necessary 
to satisfy the boundary conditions. 
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The type of solution, can also tie altered using only certain types 
of terns. For instance, if it is desired to obtain a solution for -which 
the flow leaves the body tangent to the surface, at the edge of the "base, 
such as Is usually the case in a viscous flow for todies at small angles of 
attach, rather than a solution for which the flow tri.ll negotiate the sharp 
"base edge and continue down the "base to a rear stagnation point, then only 
those terms for which the vorticlty goes to zero at the tase edge should he 
used. Fording the flow to leave the "body surface tangentially is similar 
to the Kutta condition used to define the amount of circulation around an 
airfoil. The reasons fqr such a constraint are purely empirical in nature 
and will only represent the physical nature of the flow at certain angles 
of attach. At large angles of attack the flow would he expected to leave the 
surface on the leeward side at some unknown forward longitudinal station. 

The strength of the discrete vortex ^ used to represent the vorticlty 
over an incremental distance Ji along a meridian line on the fanpod 
surface is equal to the integral of the vorticlty over the distance ^ • . 

Then /J can he written in matrix notation from equation (l); 


<t> ][W>} 


(2) 


where 


M = 






• Oj GTC. 


* . 


•j 


and is the number of discrete vortices on the fanpod in the chordwise 

direction. 

It should he noted that there is one equation; such as equation (2), 
for each roll angle. The discrete horseshoe type vortex strengths /7 can 
he transformed into a set of equivalent quadrilateral vortex strengths 
hy means of the following expression. 

M = [ K i ]M 
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■where 
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The relationship between the horseshoe nnrl quadrilat e ral vortex 
systems is best illustrated hy the following schematic of the vortex grid 



The quadrilateral vortex strengths ore then written as a function of the 
unknown vorticity function coefficients F (<©) hy substituting equation ( 2 ) 
into equation (3). 


fa*'} - *° [ K ]["^{][ ? ]fa’] 
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Since there is one equation, such as equation (U), for each roll angle ® 
the matrix equation for o sr & < 36O is given by; 
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where 
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APPENDIX I 


CROSS PLOW DUS TO LIFT 


The cross flow due to lift is a combination of that due to the ehordwise 
component of bound vorticity along constant percent chord lines and that due 
to the trailing vorticity. The spanwise component of hound vorticity is 
given hy equation (50) of Appendix p. 


. _ 4 b ^ f. 

~ C C c 


£cozr<s>/-z ■+ 


r 

A»fe/ 








/jsy> + &/) y 

Jsy^>fOs>- tyJf 




CD 


The trailing vorticity is given by; 
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therefore; the cross flow due to lift is then given by; 
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•where / is the local sweep of the constant percent chord line. 

Hhe span-wise functions 'f , & , % , and < for the wing are defined - 
in Appendix F. Their derivatives are given "below. 
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The derivatives of the ) functions for breaks in the leading 

or trailing edges are obtained from the following relations. 

For c> 


' ^$^<0 X (8) 
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Also, the derivative of f 0 r partial span leading and trailing 

edge flaps. 


^ ^ . **z.- > 

( 12 ) 


■where , %• , and are the span stations where the break in the 

leading or trailing edge occurs, where the partied span flap "begins inboard, 
and where the partial span flap ends outboard, respectively. The 
derivatives ^ and used in equations (8) through (IS) 

are derived below. The variable ' is a dummy variable. Since mCG tp*!? 
defined in equation (33) of ^Appendix A, is given as; 
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or 
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The expression for is given as equation (^3) of Appendix A. 
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Due to the singular nature of equation (21 ) special attention must be 
given to the case vhen ^**<?*. Eeturning to equation (20) and determining 
in the limit as < 2 ? — =>- 
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She integration with respect to <S> in equation (3) is completed tis follows; 
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and similarly; 
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therefore; eqiiation (3) becomes 
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vhere ^c** - . and *£< f. are the values of at the leading 

and trailing edgej respectively. 
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The equation for 0 was obtained from reference 6 . In 

reference 6 the effect of leading and trailing edge breaks on the local 
isobars is given by the following expression for the local sweep of the 
isobars . 




(36) 


where 
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The subscript ^ indicates values at a break in either the leading or 
trailing edge. Also, and s4> are the sweeps of the constant percent 

chord lines on the port and starboard sides of the break, respectively. The 
angle is the value of sweep to which the isobar line converges as the 

effect of the break station dies in the spanwise direction. If a span 
station is influenced by more than one break station the effect 'is obtained 
by interpolation. 
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APPENDIX J 


INDUCED DRAG 


The vortex drag arises from the momentum transported away from an 
arbitrary lifting system by the trailing vortices and is given (25) by the 
integral equation 


CO oO 

^ ~ 4 ^ y (T?*!? t -h ) dy Jz C 1 ) 

■« -oG> 

This result is usually simplified by allowing the control surface to recede 
indefinitely far behind the object. That is, a so called Trefftz plane 
analysis is performed. * 

For this condition 5^ - o and 
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For a plane potential wake, this latter result can be transformed (25) and 
(26) by the use of Green’s theorem to 
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This result is cast here into the following form by use of a change of 
variable 4^4, for numerical purposes 
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■with a numerically given function -f(f), $(?) continuous on the interval 
Co /) , and f 'C<=>) - $'u)- o have teen studied by Ementon (27) and (2o) 
using Fournier series analyses in conjunction with the condition that the . 
integral I be a minimum in order to eliminate numerical induced inaccuracies. 

The result for the current problem is 
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APPENDIX K 


SKIN FRICTION DRAG 


Several well established seraiempirical techniques for the evaluation 
of laminar and turbulent flat plate skin friction at incompressible and 
compressible speeds are used to develop a subroutine to estimate the viscous 
drag of advanced aircraft using a component buildup approach. Specified 
or flat plate natural transition point calculation options are provided for 
in conjunction with a matching of the momentum thickness to link the two 
boundary layer states. For the turbulent condition, the increase in drag 
due to distributed surface roughness is treated using uniformly distributed 
sand grain results. Component thickness effects are approximated using 
experimental data correlations for two-dimensional airfoil sections and 
bodies of revolution. 

Considerations such as separation, component interference, and discrete 
protuberances (e.g. antennas, drains, aft facing steps, etc.) must be 
accounted for separately if present or judged to be significant. 

In the following, a discussion is presented for a single component 
evaluation in order to simplify writing of the equations and eliminate 
multiple subscripting. The total result is obtained by a surface area 
weighted summation of the various component analyses. 


lamnar/transition 

A specified and flat plate natural transition option are provided in 
the program. The principal function of the calculation is to provide the 
conditions required to initialize the turbulent solution. In particular, 
the transition point length and momentum thickness Reynolds numbers are 
required . 
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This solution is based on the laminar Blasius result (29, chapter VII ) 
in conjunction with Eckert’s (30) compressibility transformation. This option 
permits an assessment of the reduction in skin friction drag if laminar flow 
can be maintained for the specified extent. It does not establish the likeli- 
hood that such a condition will be realized in practice or to what extent 
with the following exception. 

FIAT PLATE NATURAL 
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These techniques have been taken from Schlicting (29, chapters XVI and 
XVII ). The functions of fq and Fq which relate the effect of freestream 
turbulence intensity as the difference between the transition and critical 
Reynolds number and the effect of roughness on smooth flat plate transition 
are figures 16.21 and a flat plate normalized version of 17 * 28 , respectively. 

The flat plate natural transition option is an oversimplified description 
of what actually takes place as a result of the neglect of pressure gradient 
and three dimensional crossflow vorticity effects. 


TURBULENT 


A smooth and distributed rough surface option have been provided in the 
analysis. In either case, the solution is initialized by matching the 
momentum thickness at the transition point produced by the laminar/transition 
solution. That is, an effective origin (commonly referred to as a virtual 
origin) is established for the turbulent analysis. 
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where 


~ ^*77ZS!sO ( 2 °) 

' - C -2 ^-a)/^ r4^ ( 21 ) 
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The compressible turbulent flat plate method used here is that proposed 
by Van Driest (3l) based on the Von Karman mixing length hypothesis in . con- 
junction with the Squire-Young formulation for profile drag (29, chapter XXIV) 
as applied to a flat plate. 
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The turbulent flat plate method used here is that of Schllcting (29 ? 
chapter XXI ) which is based on a transposition of Hikuradse's densely packed 
sand grain roughened pipe data* The effect of compressibility is due to the 
reduction in density at the wall as proposed by Goddard (32). The selection 
of the equivalent sand grain roughness for a given manufacturing surface 
finish is made with aid of Table I which was taken from Clutter (3*0. 
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Type of .Surface 


TABLE I 


Equivalent Sand Roughness 
k (inches) 


Aerodynamically smooth 

Polished metal or wood 

Natural sheet metal 

Smooth matte paint, carefully applied 

Standard camouflage paint, average application 

Camouflage paint, mass -production spray 

Dip-galvanized metal surface . 

Natural surface of cast iron 


0 

0,02 - 0.08 x io -3 
0.16 x 10" 3 
0.25 x 10"3 
0 .1*0 x lo -3 
1.20 x 10“ 3 
6 X 10-3 
10 x 10“ 3 


THICKNESS CORRECTIONS 


The foregoing evaluations produce an estimate of the shearing forces on 
a flat plate (at zero angle of attack) for a variety of conditions. As an 
actual aircraft has a nonvahishing thickness, an estimate of pressure gradient 
effects on skin friction and boundary layer displacement pressure drag losses 
is required. A common procedure for accomplishing this and the one which 
will he used here is based on nonlifting experimental correlations for 
symmetric two dimensional airfoils and axisyrametric bodies. The following 
relations derived by Hoerner ( 33 ? chapter VI ) are used respectively. 


* ^ 
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Hoerner recommends K]j.“ 2 for airfoils with maximum thickness at 
30 percent chord and K]_ = 1.2 for NACA 61 * and 65 series airfoils. In this 
regard, the best information available to an analyst for his particular 
contour should be used. This is especially true for modern high performance 
shapes such as the supercritical airfoil, etc. 


2 1*1 



TOTAL VISCOUS DRAG 


An aircraft total viscous drag coefficient is estimated by a sum of the 
preceding analysis over all components (i.e. wing, fuselage, vertical tail, 
etc.). That is 

* 

^ Cp -> ^ ^ 

The component length used in the calculation of the skin friction 
coefficient is the mean aerodynamic chord for planar components or segments 
thereof and the physical length for bodies and nacelles. 
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recovery factor 

unit Reynolds number, Uco/V^, ft ~ 1 
Reynolds number based on [ ] 
gas constant, 1716 ft^ °R /sec 
surface area - ft- 

4 

static temperature - °R 
recovery temperature ratio, 1+ r Til Mco 
airfoil thickness ratio 
velocity, s/^RT 1 M - ft/sec 

longitudinal distance from beginning of component * ft 

ratio of Specific heats 

boundary layer momentum thickness , ft 
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p 

p 

r 

i 

Inst 

j 

i 

r 

TRAN 
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absolute viscosity, lbf-sec/ft 3 
kinematic viscosity, P/p ft^/sec 
density, P/iftT ,Ibf- sec 2 /ft^ 
turbulence intensity 


SUBSCRIPTS 
iteration number 
instability 
component rfumber 
recovery 

transition point 
free stream 


SUPERSCRIPTS 


property based on effective origin 
Eckert reference temperature condition 
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